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Abstract 

Gravitational-wave experiments with interferometers and with resonant masses can 
search for stochastic backgrounds of gravitational waves of cosmological origin. We 
review both experimental and theoretical aspects of the search for these backgrounds. 
We give a pedagogical derivation of the various relations that characterize the response 
of a detector to a stochastic background. We discuss the sensitivities of the large interfer- 
ometers under constructions (LIGO, VIRGO, GEO600, TAMA300, AIGO) or planned 
(Avdanced LIGO, LISA) and of the presently operating resonant bars, and we give 
the sensitivities for various two-detectors correlations. We examine the existing limits 
on the energy density in gravitational waves from nucleosynthesis, COBE and pulsars, 
and their effects on theoretical predictions. We discuss general theoretical principles 
for order-of-magnitude estimates of cosmological production mechanisms, and then we 
turn to specific theoretical predictions from inflation, string cosmology, phase transi- 
tions, cosmic strings and other mechanisms. We finally compare with the stochastic 
backgrounds of astrophysical origin. 
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1 Introduction 



A large effort is presently under way in the scientific community for the construction of 
the first generation of large scale interferometers for gravitation wave (GW) detection. The 
two LIGO detectors are under construction in the US, and the Italian-French interferometer 
VIRGO is presently under construction near Pisa, Italy. The detectors are expected to start 
taking data around the year 2002, with a sensitivity where we might find signals from known 
astrophysical sources. These detectors are also expected to evolve into second generation 
experiments, with might better sensitivities. At the same time, two other somewhat smaller 
interferometers arc being bTult: GEO600 in Germany and TAMA300 in Japan, and they will 
be important also for developing the techniques needed by second generation experiments. 
To cover the low-frequency region, where many interesting signals are expected, but which 
on the Earth is unaccessible because of seismic noise, it is planned to send an interferometer 
into space: this extraordinary experiment is LISA, which has been selected by the European 
Space Agency has a cornerstone mission in his future science program Horizons 2000 Plus, 
and might fiight somewhere around 2010-2020, with a sensitivity levels where many sources 
are virtually guaranteed. Finally, although at a lower sensitivity level, cryogenic resonant 
bars are in operation since 1990, and more sensitive resonant masses, e.g. of spherical shape, 
are under study. 

A possible target of these experiments is a stochastic background of GWs of cosmological 
origin. The detection of such a background would have a profound impact on early Universe 
cosmology and on high-energy physics, opening up a new window and exploring very early 
times in the evolution of the Universe, and correspondingly high energies, that will never be 
accessible by other means. The basic reason why relic GWs carry such a unique information 
is that particles which decoupled from the primordial plasma at time t ~ tdec, when the 
Universe had a temperature Tdec, give a 'snapshot' of the state of the Universe at T '-^ T^ec- 
All informations on the Universe when the particle was still in thermal equilibrium has 
instead been obliterated by the successive interactions. The weaker the interaction of a 
particle, the higher is the energy scale when they drop out of thermal equilibrium. Hence, 
GWs can probe deeper into the very early Universe. 

To be more quantitative, the condition for thermal equilibrium is that the rate T of the 
processes that mantain equilibrium be larger than the rate of expansion of the Universe, as 
measured by the Hubble parameter H. The rate (in units ^ = c = 1) is given by F = na\v\ 
where n is the number density of the particle in question, and for massless or light particles 
in equilibrium at a temperature T, n ~ l^l ~ 1 is the typical velocity and a is the cross- 
section of the process. Consider for instance the weakly interacting neutrinos. In this case 
the equilibrium is mantained, e.g., by electron-neutrino scattering, and at energies below 
the W mass a ~ G^{E^) ~ G'^T^ where Gp is the Fermi constant and {E'^) is the average 
energy squared. The Hubble parameter during the radiation dominated era is related to 
the temperature by H ^ T^/Mpi. Therefore 



Even the weakly interacting neutrinos, therefore, cannot carry informations on the state 
of the Universe at temperatures larger than approximately 1 MeV. If we repeat the above 




(1) 
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computation for gravitons, the Fermi constant Gp is replaced by Newton constant G = 
1/Mpj, where Mpi ~ 10^^ GeV is the Planck mass. At energies below Mpi, 



The gravitons are therefore decoupled below the Planck scale. It follows that relic gravita- 
tional waves are a potential source of informations on very high-energy physics. Gravita- 
tional waves produced in the very early Universe have not lost memory of the conditions 
in which they have been produced, as it happened to all other particles, but still retain 
in their spectrum, typical frequency and intensity, important informations on the state of 
the very early Universe, and therefore on physics at correspondingly high energies, which 
cannot be accessed experimentally in any other way. It is also clear that the property of 
gravitational waves that makes them so interesting, i.e. their extremely small cross section, 
is also responsible for the difficulties of the experimental detection. 

In this review we will examine a number of experimental and theoretical aspects of the 
search for a stochastic background of GWs. The paper is organized as follows. We start 
from the experimental side; in sects. 2 to 4 we present the basic concepts for characterizing 
a stochastic background and for describing the detector. This seems to be a subject where 
everybody has its own definitions and notations, and we made an effort to clean up many 
formulas appearing in the literature, to present them in what we consider the clearest form, 
and to give the most straightforward derivation of various relations. 

In sect. 5 we examine the various detectors under construction or already operating; 
we show the sensitivity curves of LIGO, VIRGO, GEO600, Advanced LIGO, LISA and 
NAUTILUS (kindly provided by the various collaborations), and we discuss the limiting 
sources of noise and the perspectives for improvement in the near and mid-term future. 

To detect a stochastic background with Earth-based detectors, it turns out that it is 
mandatory to correlate two or more detectors. The sensitivities that can be reached with 
two-detectors correlations are discussed in sect. 6. We will show graphs for the overlap 
reduction functions between VIRGO and the other major detectors, and we will give the 
estimate of the minimum detectable value of the signal. 

In sect. 7 we examine the existing bounds on the energy density in GWs coming from 
nucleosynthesis, COBE, msec pulsars, binary pulsar orbits, pulsar arrays. 

Starting from sect. 8 we move toward more theoretical issues. First of all, we discuss gen- 
eral principles for order of magnitude estimate of the characteristic frequency and intensity 
of the stochastic background produced by cosmological mechanisms. With the very limited 
experimental informations that we have on the very high energy region, MGUT~^^^Pb 
and on the history of the very early Universe, it is unlikely that theorists will be able to 
foresee all the interesting sources of relic stochastic backgrounds, let alone to compute their 
spectra. This is particularly clear in the Planckian or string theory domain where, even if 
we succeed in predicting some interesting physical effects, in general we cannot compute 
them reliably. So, despite the large efforts that have been devoted to understanding possible 
sources, it is still quite possible that, if a relic background of gravitational waves will be 
detected, its physical origin will be a surprise. In this case it is useful to understand what 
features of a theoretical computation have some general validity, and what are specific to 
a given model, separating for instance kinematical from dynamical effects. The results of 




(2) 
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this section provide a sort of benchmark, against which we can compare the results from 
the specific models discussed in sect. 9 and 10. 

In sect. 9 we discuss the most general cosmological mechanism for GW production, 
namely the amplification of vacuum fluctuations. We examine both the case of standard 
inflation, and a more recent model, derived from string theory, known as pre-big-bang 
cosmology. In sect. Jl^ we examine a number of other production mechanisms. 

Finally, in sect. |l^ we discuss the stochastic background due to many unresolved astro- 
physical sources, which from the point of view of the cosmological background is a 'noise' 
that can compete with the signal (although, of course, it is very interesting in its own right). 



2 Characterization of stochastic backgrounds of GWs 

A stochastic background of GWs of cosmological origin is expected to be isotropic, station- 
ary and unpolarized. Its main property will be therefore its frequency spectrum. There 
are different useful characterization of the spectrum: (1) in terms of a (normalized) energy 
density per unit logarithmic interval of frequency, /igr2gw(/); (2) in terms of the spectral 
density of the ensemble avarage of the Fourier component of the metric, Sh{f)', (3) or in 
terms of a characteristic amplitude of the stochastic background, hc{f). In this section we 
examine the relation between these quantities, that constitute the most common description 
used by the theorist. 

On the other hand, the experimentalist expresses the sensitivity of the apparatus in 
terms of a strain sensitivity with dimension Hz"^/^, or thinks in terms of the dimensionless 
amplitude for the GW, which includes some form of binning. The relation between these 
quantities and the variables Sh{f), hQ^g„{f), hdf) will be discussed in sect. 3. 



2.1 The energy density Qg„{f) 



The intensity of a stochastic background of gravitational waves (GWs) can be characterized 
by the dimensionless quantity 

Pc 



^^gw(/) 



(3) 



dlogf 

where pgw is the energy density of the stochastic background of gravitational waves, / is 
the frequency (w = 27r/) and pc is the present value of the critical energy density for closing 
the Universe. In terms of the present value of the Hubble constant Hq, the critical density 
is given by 

(4) 



Pc 



8ttG 



The value of Hq is usually written as Hq = Hq x 100 km/(sec-Mpc), where Hq parametrizes 



the existing experimental uncertainty. Ref. |20£] gives a value 0.5 < Hq < 0.85. In the last 
few years there has been a constant trend toward lower values of Hq and typical estimates 
are now in the range 0.55 < ho < 0.60 or, more conservatively, 0.50 < ho < 0.65. For 
instance ref. |226], using the method of type lA supernovae, gives two slightly different 
estimates ho = 0.56 ± 0.04 and ho = 0.58 ± 0.04. Ref. [247], with the same method, finds 
Hq = 0.60 lb 0.05 and ref. |165|, using a gravitational lens, finds Hq = 0.51 it 0.14. The spread 



of values obtained gives an idea of the systematic errors involved. 
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It is not very convenient to normalize to a quantity, pc, which is uncertain: this 
uncertainty would appear in all the subsequent formulas, although it has nothing to do 
with the uncertainties on the GW background itself. Therefore, we rather characterize the 
stochastic GW background with the quantity /ig$7gw(/), which is independent of h^. All the- 
oretical computations of primordial GW spectra are actually computations of dpg^/dlog f 
and are independent of the uncertainty on Hq. Therefore the result of these computations 
is expressed in terms of /igf^gw, rather than of f^gwQ 

2.2 The spectral density Shif) and the characteristic amphtudes hc{f) 

To understand the effect of the stochastic background on a detector, we need however to 
think in terms of amplitudes of GWs. A stochastic GW at a given point x = can be 
expanded, in the transverse traceless gauge, a^ 




where hA{—f,Cl) = h\{f,(l). is a unit vector representing the direction of propagation 
of the wave and dQ = d cos Odcj). The polarization tensors can be written as 

6afe(^) = marhb-hahb, 

^abi^) = rhafih + narhh, (6) 

with m, n unit vectors ortogonal to (l and to each other. With these definitions, 

e^,(A)e^>^(J7) = 25^^' . (7) 

For a stochastic background, assumed to be isotropic, unpolarized and stationary (see |^, |l^ 
for a discussion of these assumptions) the ensemble average of the Fourier amplitudes can 
be written as 

{h\{f, n)hA'if', n')) = 5{f - f')^s\n, n')6AA'lsh{f) , (8) 

where 6^{(l,(l') = 5{(p — (p')S{cos9 — cos 6'). The spectral density Sh{f) defined by the 
above equation has dimensions Hz~^ and satisfies Sh{f) = Sh{—f)- The factor 1/2 is 
conventionally inserted in the definition of Sh in order to compensate for the fact that the 
integration variable / in eq. (|5|) ranges between — oo and +00 rather than over the physical 
domain < / < 00. The factor 1/{4tt) is a choice of normalization such that 

I dCidh' {h*Aif,h)hA'{f',n')) = 6{f - f')6AA'lshif) . (9) 

^This simple point has occasionally been missed in the literature, where one can find the statement that, 
for small values of Ho, figw is larger and therefore easier to detect. Of course, it is larger only because it 
has been normalized using a smaller quantity. 

^Our convention for the Fourier transform are g(f) ~ dt exp{2Tvift}g{t), so that g{t) = 

J^^df e^p{~27vift}g{f). 
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Using eqs. (|5|,^) we get 

/oo rf=OQ 
dfShU) = A d{logf) fShif). (10) 

-oo Jf=0 

We now define the characteristic amplitude hdf) from 

/•/=oo 

Mt)h''\t)) = 2 d{logf)hl{f). (11) 

Jf=0 

Note that hdf) is dimensionless, and represents a characteristic value of the amplitude, per 
unit logarithmic interval of frequency. The factor of two on the right-hand side of eq. ( pi] ) 
is part of our definition, and is motivated by the fact that the left-hand side is made up 
of two contributions, given by {h*^hj^) and {h*^hx)- In a unpolarized background these 
contributions are equal, while the mixed term {h*^hx) vanishes, eq. (|8|). 
Comparing eqs. (|lO|) and (|Tl|), we get 

hl{f) = 2fSdf)- (12) 

We now relate hdf) and hQ^lg„{f). The starting point is the expression for the energy 
density of gravitational waves, given by the 00-component of the energy-momentum tensor. 
The energy-momentum tensor of a GW cannot be localized inside a single wavelength (see 
e.g. ref. [p01[| , sects. 20.4 and 35.7 for a careful discussion) but it can be defined with a 
spatial averaging over several wavelengths: 

Pgw = ^(/la6/l'^'). (13) 

For a stochastic background, the spatial average over a few wavelengths is the same as a 
time average at a given point, which, in Fourier space, is the ensemble average performed 
using eq. (^). We therefore insert eq. (^) into eq. (13) and use eq. (P). The result is 

Pgw = ^ J^^~^°^ d(log /) f{2^f)^Sh{f) , (14) 

SO that 

rin.. TT 

;fS,if) . (15) 



dlogf 2G 

Comparing eqs. ([T5|) and (12) we get the important relation 

dpgw 



dlogf AG 
or, dividing by the critical density pc, 



-JT^fXif), (16) 



97r2 

^^U) = ^fhlU)- (17) 



Using eq. (12) we can also write 



^^gw(/) = ^/'S,(/). (18) 
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Some of the equations that we wih find below are very naturally expressed in terms of 
hQQgviif)- This will be true in particular for all theoretical predictions. Conversely, all 
equations involving the signal-to-noise ratio and other issues related to the detection are 
much more transparent when written in terms of Sh{f)- Eq. ( p!8| ) will be our basic formula 
for moving between the two descriptions. 

Inserting the numerical value of Hq, eq. ([ITD gives 



h,if) 1.263 X 10-1^ (^^) \lhln^{f) . (19) 

Actually, hc{f) is not yet the most useful dimensionless quantity to use for the comparison 
with experiments. In fact, any experiment involves some form of binning over the frequency. 
In a total observation time T, the resolution in frequency is A/ = l/T, so one does not 
observe /io^gw(/) but rather 

d(log /) hln^^if) ^ -^hln^^if) , (20) 

and, since hQ^lg^{f) ~ h^f), it is convenient to define 

h,(f,Af) = h,{f)^-^j . (21) 

Using 1/(1 yr) ~ 3.17 x 10"*^ Hz as a reference value for A/, and 10~^ as a reference value 
for /igf^gw, eqs. (|l9|) and (|2l|) give 

or, with reference value that will be more useful for LISA, 

A/) = 7^in X 10- rJ^f C^^] ( ^) . (23) 

■'^ V / / V / V3.17 X 10-8Hzy ^ ^ 

Finally, we mention another useful formula which expresses hQQg^{f) in terms of the number 
of gravitons per cell of the phase space, n{x,k). For an isotropic stochastic background 
n{x,k) = Uf depends only on the frequency / = |fc|/(27r), and 



Therefore 



and 



^ 271 fn f = levr^ d{\og f)fnf. (24) 



dlogf 

4 



3^6 (^) (^) . (26) 
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This formula is useful in particular when one computes the production of a stochastic 
background of GW due to amplification of vacuum fluctuations, since the computation of 
the Bogoliubov coefficients (see sect. gives directly nj. 

As we will discuss below, to be observable at the LIGO/VIRGO interferometers, we 
should have at least /igfigw ~ 10~^ between 1 Hz and 1 kHz, corresponding to of order 
10^^ at 1 kHz and ~ 10^^ at 1 Hz. A detectable stochastic GW background is therefore 
exceedingly classical, 1. 



3 The response of a single detector 

3.1 Detector tensor and detector pattern functions 

The quantities Ogw(/), /ic(/)i 5'/i(/) discussed above are all equivalent characterization of 
the stochastic background of GWs, and have nothing to do with the detector used. Now 
we must make contact with what happens in a detector. The total output of the detector 
S{t) is in general of the form 

S{t) = s{t) +n{t) (27) 

where n{t) is the noise and s{t) is the contribution to the output due to the gravitational 
wave. For instance, for an interferometer with equal arms of length L in the x — y plane, 
s{t) = [5Lx{t) — 5Ly{t))/L, where SL^^y are the displacements produced by the GW. The 
relation between the scalar output s{t) and the gravitational wave habit) in the transverse- 
traceless gauge has the general form |117] , |113| 

s{t) = D''''Kb{t) . (28) 

D"-^ is known as the detector tensor. 
Using eq. (^, we get 

/oo r 
df / dn hAif, n)e-'-'^^'D-'ei{n) . (29) 

/i=+,x -°° 

It is therefore convenient to define the detector pattern functions Fa{^), 

= D^'e^.m (30) 

SO that 

/oo r 
df / dCl JiAif, n)FA{Cl)e-^^'f' , (31) 

and the Fourier transform of the signal, s(/), is 

~<f)= E fdnhAif,Cl)FAiCl). (32) 

A=+,x 

The pattern functions Fa depend on the direction = (0, (p) of arrival of the wave. Fur- 
thermore, they depend on an angle (hidden in e^^) which describes a rotation in the plane 
orthogonal to (l, i.e. in the plane spanned by the vectors m, n, see eq. (|6[). Once we have 
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made a definite choice for the vectors rh, h, we have chosen the axes with respect to which 
the + and x polarizations are defined. For an astrophysical source there can be a natural 
choice of axes, with respect to which the radiation has an especially simple form. Instead, 
for a unpolarized stochastic background, there is no privileged choice of basis, and the angle 
■i/; must cancel from the final result, as we will indeed check below. 

For a stochastic background the average of s{t) vanishes and, if we have only one de- 



tector, the best we can do is to consider the average of s (t). Using eqs. (pC|) and 



oo 1 rco 



{s'it)) = F I df\sH{f) = F I dfShif), (33) 



where 

(in 







^^J^H i^^(^^,V')i^^(^^,V'). (34) 

Note that, while the value of {s'^{t)) is obtained summing over all stochastic waves coming 
from all directions Cl, the factor l/(47r) in eq. (^) produces an average of + F^ over the 
directions. The factor F gives a measure of the loss of sensitivity due to the fact that the 
stochastic waves come from all directions, compared to the sensitivity of the detector for 
waves coming from the optimal direction. 

To compute the pattern functions, we need the explicit expressions for the vectors rh, h 
that enters the definition of eq. (|6|). We use polar coordinates centered on the detector. 
Then 

= (sin 9 sin (j), sin 9 cos (j), cos 6) (35) 
and a possible choice for m, n is 

h = {— cos9 sincp, — cos cos (I),sm9) , rfi = (— cos iji, sin (/>, 0) . (36) 

The most general choice is obtained with the rotation 

h ^ n cos ip + rh sin ip , 

rh — > —h sin ip + rh cos tp , (37) 

so that 

fi = (— cos sin i;^ cos — cos sin -i/), — cos ^ cos (/> cos "0 + sine/) sin -0, sin cos -0) , 
m = (— cos (f) cos -0 + cos 9 sin (f) sin tp, sin (p cos tp + cos 9 cos (p sin ip, — sin 9 sin ^) . (38) 

We now compute the explicit expression for F^ and F in the most interesting cases. 
3.1.1 Interferometers 

For an interferometer with arms along the ii and v directions (not necessarily orthogonal) 

D^^ = ^ (^u^u^ - v%'') . (39) 
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Using the definition of Fa, eq. (|30|), together with eq. (^) for e^^ and with the above 
expressions for rh, h one obtains, restricting now to an interferometer with perpendicular 



arms [|lT|, g3|, |4|, 

The factor F is then given by 



(1 + cos^ 6) cos 2(j) cos 2'0 — cos 9 sin 2^ sin 2^/^ ; 
(1 + cos^ 9) cos 2(f) sin 2'0 + cos 9 sin 2(/) cos 2tp . 



(40) 



(41) 



Note also that {s'^it)) depends on the Fa only through the combination F^ + F^ which is 
independent of the angle ip, as can be checked from the explicit expressions, in agreement 
with the argument discussed above. The same will be true for resonant bars and spheres. 
Therefore below we will often use the notation FA{i^) instead of Fa{^,iP)- 

It is interesting to see how these results are modified if the arms are not perpendicular. 
If a is the angle between the two arms, we find 



F+(^, ■0) = (sin a) 

Fx(^,</',V) = (sino) 
Therefore 



1 + cos^ 9) sin(a + 2</)) cos 2^^ + cos 9 cos(a + 24>) sin 2^^ 
1 + cos^ 9) sin(a + 2cj)) sin 2^^ — cos 9 cos(a + 20) cos 2^/> 

r = — sm a . 



(42) 



(43) 



For a = 7t/2 we recover the results (4C,41) and the sensitivity is maximized, while for a = 
the sensitivity of the interferometric detection of course vanishes. 



3.1.2 Cylindrical bars 

For a cylindrical bar with axis along the direction /, one has instead 



(44) 



Since -D°* is contracted with the traceless tensor hab, it is defined only apart from terms 
~ 5°'', so that we can equivalently use the traceless tensor 



D 



ab 



1 



:ab 



(45) 



To compute Fa it is convenient, compared to the case of the interferometer, to perform a 
redefinition ifj — > '0 + 7r/2, and define 9 as the polar angle measured from the bar direction /. 



(of course, when we will discuss interferomer-bar correlations in sect. 4.3, we will be careful 
to use the same definitions in the two cases). Then 



F+{9A,^) 



sin^ 9 cos 2^lJ , 
sin^ 9 sin 2^/^ . 



(46) 



11 



and the angular factor F is 



(47) 



A=+,x 



3.1.3 Spherical detectors 



Finally, it is interesting to give also the result for a detector of spherical geometry. For a 
spherical resonant-mass detector there are five detection channels |157] ] corresponding to 
the five degenerates quadrupole modes. The functions D""^ and Fa for each of these channel 



are given in ref. [261|: one defines the real quadrupole spherical harmonics as 



I20 
1 



V2 



(y2,-i-i2,+i) 



^1. 

Y2s 



^{Y2,-l+Y2,+l) 
^{Y2,-2+Y2,+2) 



(48) 



72 



{Y2,-2 - Y2,+2) ■ 



The pattern functions associated to the 0, Ic, Is, 2c, 2s channels are 



F+,2c{9,(t>'. 



■ sm 



— sin 29 cos ( 
2 

-— sin 20 sin < 
2 



^(1 + cos^6')cos20 
^(1 + cos2 6')sin20 



Fx,o{e,<P) -- 
Fx,ic{0,<P) 
i^x,is(e,0) 

Fx,2c{9^<P) 

FxM^^^) 







sm f sm (p 
sin 9 cos (j) 
cos 9 sin 2(j) 
cos 9 cos 20 . 



(49) 



Here the angle ip has been set to zero; as discussed above, this corresponds to a choice of the 
axes with respect to which the polarization states are defined. With this choice, at 9 = 0, 
only Fa,2s,Fa,2c are non- vanishing, consistently with the fact GWs have elicities ±2. The 
same value as for interferometers, F = 2/5, is obtained for the sphere in the channel m = 0, 
as well as for the channels m = ±1 and m = ±2. 

3.2 The strain sensitivity hj 

The ensemble average of the Fourier components of the noise satisfies 



{n*imf')) = 6{f-f')lsM)- 



(50) 
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The above equation defines the functions Sn{f), with /) = Snif) and dimensions 
Hz^^. The factor 1/2 is again conventionally inserted in the definition so that the total 
noise power is obtained integrating Sn{f) over the physical range < / < cxd, rather than 
from — oo to oo, 

{n\t))= / dfSM)- (51) 







The function Sn is known as the square spectral noise density^ 

Equivalently, the noise level of the detector is measured by the strain sensitivity h f 

hf^^/sM), (52) 

where now / > 0. Note that ^/ is linear in the noise and has dimensions Hz~^/^. 

Comparing eqs. ( |5ll ) and eq. (p3|), we see that in a single detector a stochastic background 
will manifesting itself as an eccess noise, and will be observable at a frequency / if 

Sh(.f)>jSM)- (53) 

where F is the angular efficency factor, eq. (34). Using eq. ( [I8| ) and S'„(/) = h'j, we can 
express this result in terms of the minimum detectable value of /igilgw, as 

^W(/)4-10-(^)'(^5.3^)'. (54) 

In sect. 4 we will use this result to compute the sensitivity of various single detectors to a 
stochastic background. 



4 Two-detectors correlation 

4.1 The overlap reduction functions 7(/),r(/) 

To detect a stochastic GW background the optimal strategy consists in performing a cor- 
relation between two (or more) detectors, since, as we will see, the signal is expected to be 
far too low to exceed the noise level in any existing or planned single detector (with the 
exception of the planned space interferometer LISA) . 

The strategy for correlating two detectors has been discussed in refs. [^00| , 74, 113, 257 , 



^. We write the output Si{t) of the i— th detector as Si{t) = Si{t) + ni{t), where i = 1,2 
labels the detector, and we have to face the situation in which the GW signal Si is much 
smaller than the noise n^. At a generic point x we rewrite eq. (^) as 

habit, x)= r df [dClhA{f,Cl)e-^'''^^'-^ '^/'^^ei{n), (55) 



4 



Unfortunately there is not much agreement abou t not ations in the hterature. The square spectral noise 



density, that we denote by Sn{f) following e.g. ref. 1113], is called P{f) in ref. |8|. Other authors use the 



notation Sh{f ), which we instead reserve for the spectral density of the signal. To make things worse, S„ is 
sometime defined with or without the factor 1/2 in eq. (pOl). 
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so that 



/CO P 
-oo J 



where F/^ are the pattern functions of the i-th detector. The Fourier transform of the scalar 
signal in the i-th detector, Si{f), is then related to hA{f,Cl) by 

^^if)= E fdnhA{f,n)e'^'^^-^'/'F,^{Cl). (57) 
A=+,x 

We then correlate the two outputs defining 

rT/2 fT/2 

Si2= dt dt' Si{t)S2{t')Q{t-t'), (58) 

J-T/2 

where T is the total integration time (e.g. one year) and Q a real filter function. The 
simplest choice would be Q{t — t') = 5{t — t'), while the optimal choice will be discussed in 
sect. [4.2[ At any rate, Q{t — t') falls rapidly to zero for large \t — t'\; then, taking the limit 
of large T, eq. (|58|) gives 

/ + 00 
df Sl{f)S2{f)Q{f) . (59) 
-oo 

Taking the ensemble average, the contribution to {S12) from the GW background is 



{su) ^ / df{sUf)s2{f))Q{f)= df dCidn'e^^''''^-''-^'^ 



X 

xJ2Ff{n)F^'in'){h*Aif,n)hA'{f,n'))Q{f) = 

A,A' 

df\Sh{f)T{f)Q{f), (60) 
where we have used eq. (^ and (5(0) = /Zx'y2 dt = T. In the last line we have defined 



exp|27ri/f7 - — I , (61) 



Y^F,^{n)F,^{n) 

. A 

where Ax is the separation between the two detectors. Note that J2a ^1^2 is independent 
of ip even if the two detectors are of different type, e.g. an interferometer and a bar, as can 
be checked from the explicit expressions, and therefore again we have not written explicitly 
the il) dependence. 
We introduce also 

dil 



'12 



/ ^E^i'(^)^2^(^)laligned, (62) 

where the subscript means that we must compute F12 taking the two detectors to be per- 
fectly aligned, rather than with their actual orientation. Of course if the two detectors are 
of the same type, e.g. two interferometers or two cylindrical bars, F12 is the same as the 
constant F defined in eq. (f 
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The overlap reduction function 7(/) is defined by fTl, 113 



T(f) 

lU) = -j^ ■ (63) 

This normahzation is useful in the case of two interferometers, since F12 = 2/5 aheady 
takes into account the reduction in sensitivity due to the angular pattern, already present 
in the case of one interferometer, and therefore 7(/) separately takes into account the effect 
of the separation Ax between the interferometers, and of their relative orientation. With 
this definition, 7(/) = 1 if the separation Arc = and if the detectors are perfectly aligned. 

This normalization is instead impossible when one considers the correlation between 
an interferometer and a resonant sphere, since in this case for some modes of the sphere 
F12 = 0, as we will see below. Then one simply uses r(/), which is the quantity that enters 
directly eq. (|60|). Furthermore, the use of r(/) is more convenient when we want to write 
equations that hold independently of what detectors (interferometers, bars, or spheres) are 
used in the correlation. 



4.2 Optimal filtering 



As we have seen above, in order to correlate the output of two detectors we consider the 
combination 

'S'12 = 



df Slif)S2{f)Q{f) 



(64) 



We now find the optimal choise of the filter function Q{f), [^OOl , 74, 113, |257| , |8| following 
the discussion given in ref. Since the spectral densities are defined, in the unphysical 
region / < 0, by Sh{f) = Sh{-f),Sn{f) = 5'„(-/), we also require Q{f) = Q{-f) and 
Qif) real (so that Q(/) = Q*{-f) and Q{t) is real). 

We consider the variation of S'12 from its average value. 



N = S- 



12 



{S- 



12 



(65) 



By definition (N) = 0, while 



{s!2) - {S12)' 



dfdf'Q{f)Q*{f) X 



{si{f)S2{f)Si{f')s*2{n) - {snf)S2{f)){s*2{f)Si{f)) 



(66) 



We are now interested in the case rij ^ hi. Then Si — rif, if the noise in a single detector 
has a Gaussian distribution and the noise in the two detectors are uncorrelated, then this 
becomes 

(67) 
(68) 



(iV2) / dfdf'Q{f)Q*{f'){Sl{f)Si{f')){S*,{f')S2{f)) 



and 



{SKf)Si{f')) ^ {n*{m{f')) = S{f - f')lsii\f) . 



Using 5(0) = T, one therefore gets 



T 
1 



df\Qif)\'Sl{f) 
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where we have defined 

Sn{f) = HHf)SPif)Y^' ■ (70) 



We now define the signal-to-noise ratio (SNR), 



SNR 



(iV2)l/2 



(71) 



Note that we have taken an overall square root in the definition of the SNR since S12 is 
quadratic in the signal, and so the SNR is linear in the signal (this differs from the convention 
used in ref. ||8|, where it is quadratic); (S12) is the contribution to (/S'12) from the GW, given 
in eq. (|60|) . 

We now look for the function Q{f) that maximizes the SNR. For two arbitrary complex 
functions A(f), B{f) one defines the (positive definite) scalar product 

/oo 
dfA*{f)B{f)Sl{f). (72) 
-00 



Then 



and, from eq. (BH), 



Then we have to maximize 



<iV') = f(Q,Q) (73) 



M = ?(Q.f). (74) 



(SNK)«.^.x(g,^)^. ,7.) 



The solution of this variational problem is standard, 

QU) = (76) 

with c an arbitrary normalization constant. Note that, since we used r(/) rather then 7(/), 
this formula holds true independently of the detectors used in the correlations. 

It is also important to observe that the optimal filter depends on the signal that we 
are looking for, since Sh{f) enters eq. (^). This means that one should perform the data 
analysis considering a set of possible filters, chosen either in order to encompass a range of 
typical behaviours, e.g. chosing Sh{f) ~ for different values of the parameter a, or using 
the theoretical predictions discussed in later chapters. Of course this point is not relevant 
for correlations involving resonant bars, because of their narrow bandwidth. 

4.3 The SNR for two-detectors correlations 

We can now write down the value of the SNR obtained with the optimal filter: 

{Sm' = T(^,^) , (77) 
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or 



SNR: 



slif) 



nl/4 



In particular, for two interferometers r(/) = {2/5)j{f) and 



(SNR)i,tf- 



intf 



1/4 



(78) 



(79) 



For two cylindrical bars, r(/) = (8/15)7(/) and 

'l28 



(SNR)bar —bar 



.1/4 



(80) 



For the correlation between an interferometer and a cylindrical bar, again r(/) = (2/5)7(/), 
so that we get again eq. ([79[). 

In these formulas, the effect due to the fact that the two detectors will have in general 
a different orientation is fully taken into account into j{f). If however one wants to have a 
quick order of magnitude estimate of the effect due to the misalignement without performing 
each time numerically the integral that defines 7(/), one can note, from eq. (|6l|), that in 
the limit 27r fd <^ 1, where d = \Ax\ is the distance between the detectors, 



r(/) 



dn 



(81) 



where of course the pattern functions F"^ relative to the two detectors are computed with 
their actual orientation. Therefore in the limit 271 fd — > the ratio between the value of F 
that takes into account the actual orientation and the value for aligned detectors is equal to 
the ratio between F12 computed with the actual orientation and F\2 for aligned detectors. 

In the definition of 7 enters F12 computed for aligned detectors, eq. (62). Then a 
simple estimate of the effect due to the misalignement is obtained substituting in eq. ( [78|) 
F(/) ~ Fi2'y{f) where F12 is now computed for the actual orientation of the detectors and 
7(/) is the overlap reduction function for the perfectly aligned detectors. 

We apply this procedure to estimate the effect of misalignement on the SNR, in the 
correlation between an interferometer and a cylindrical bar. In the limit Infd ^ 1 we can 
also neglect the effect of the Earth curvature. Then, in the frame where the interferometer 
arms are along the x, y axes, the bar also lies in the (x, y) plane and its direction is in 
general given by / = (sin a, cos a, 0); the angle a measures the misalignement between the 
interferometer y axis and the bar. In this frame, computing the detector pattern function 
of the bar we get 



(bar) 

f 

(bar) 



COS^ 9 COS^( 



cos^ 9 cos^( 



q) cos 2-0 + sin^ 
a) sin 21^ + sin^ 



- a) cos 2-0 + cos 9 sin 2( 
q) sin 2-0 — cos 9 sin 2(( 



- q) sin 20 
a) cos 2iIj 



(82) 



To compare with eq. ( |46|) note that here 6 is measured from the z axis, and the bar is 
in the {x, y) plane, while in eq. (^6|) 6 is measured for the direction of the bar. For the 
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interferometer the functions Fa are given in eq. (pOl). A straightforward computation shows 
that 

1 



(1 + COS 6) COS 2(f) sin {(f> — a) — cos 6 cos (0 — a) 



— cos^ sin 2{(j) — a) sin 2(/) . 
Note that again the dependence on the angle ip has cancelled. Then 



'12 



47r 



COS 2a . 



(83) 



(84) 



The overall sign of F12 is irrelevant since r(/) enters quadratically in the SNR. Therefore 
we get 

SNR(a) ~ I cos 2a|i/2sjsjj^(o) . (85) 

Of course, the SNR is maximum if the bar is aligned along the y axis {ip = 0) of along the 
X axis (93 = 7r/2); the correlation becomes totally ineffective when a = 7r/4. 

Finally, it is interesting to consider the correlation between an interferometer and a 
sphere [^]. Using eq. (^9[), and eq. ( |40| ) with ip = 0, one finds in particular that for the 
correlation with the m = channel of the sphere, 



(sph)^(intf) 



A 



\/3 

— sin2 6'(l + cos^l 



I cos 20 , 



(86) 



and therefore F12 vanishes because of the integral over cp. Similarly for the channel Is, Ic 
one gets an integral over cos0cos2(/), which again vanishes. In this case of course 7(/) 
cannot be defined and one simply uses r(/). Instead, for each of the two modes 2s, 2c, 
F12 is just one half of the value for the interferometer-interferometer correlation, so that 
summing over all modes of the sphere one gets F12 = 2/5. 

4.4 The characteristic noise hn{f) of correlated detectors 



We have seen in section |2.2| that there is a very natural definition of the characteristic am- 
plitude of the signal, given by hc{f), which contains all the informations on the physical 
effects, and is independent of the apparatus. We will now associate to hc{f) a correspond- 
ing noise amplitude hn{f), that embodies all the informations on the apparatus, defining 
hc{f)/hn{f) in terms of the optimal SNR. 

If, in the integral giving the optimal SNR, eq. ([78|), we consider only a range of frequen- 
cies A/ such that the integrand is approximately constant, we can write 



SNR- 



2TA/ 



1/4 



1/4 



(87) 



where we have used eq. (|T2|). The right-hand side of eq. ( pTj ) is proportional to hc{f), and 
we can therefore define hn(f ) equating the right-hand side of eq. (|^) to hc{f)/hn{f), so 
that 



{\TMy/^ V r(/) 
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In the case of two interferometers approximately at the same site, 



r(/) 



L A 



and we recover eq. (66) of ref. | 244 |, 

hnU) = 



(4TA/)i/4 



fSnU) 



F 



1/2 



(89) 



(90) 



For comparison, note that the quantity denoted Sh{f) in ref. [244] is cahed here Sn{f)/2, 
andF = (F2) + (F2) = 2(F2). 

Prom the derivation of eq. ( p^ ) we can better understand the hmitations impUcit in the 
use of hn{f)- It gives a measure of the noise level only under the approximation that leads 
from eq. ([78|), which is exact (at least in the limit hi <^ rii), to eq. (^). This means that A/ 
must be small enough compared to the scale on which the integrand in eq. ( [78|) changes, 
so that T{f)Sh{f)/ Sn{f) must be approximately constant. In a large bandwidth this is 
non trivial, and of course depends also on the form of the signal; for instance, if /igJ7g^(/) 
is flat, as in many examples that we will find in later chapters, then Sh{f) ~ 1//^- For 
accurate estimates of the SNR at a wideband detector there is no substitute for a numerical 
integration of eq. ([78[). However, for order of magnitude estimates, eq. (p^) for hc{f) together 
eq. ( p^ ) for hn{f) are simpler to use, and they have the advantage of clearly separating the 
physical effect, which is described by hdf), from the properties of the detectors, that enter 
only in /i„(/). 

Eq. (^8|) also shows very clearly the advantage of correlating two detectors compared 
with the use of a single detector. With a single detector, the minimum observable signal, at 
SNR=1, is given by the condition Sh{f) > Sn{f)/F. This means, from eq. (|T^), a minimum 
detectable value for hdf) given by 



h 



Id 



(91) 



The superscript Id reminds that this quantity refers to a single detector. Prom eq. (| 
we find for the minimum detectable value with two interferometers in coincidence, h 
hn{f)- For close detectors of the same type, r(/) ~ F, so that 



2d 



^min(/) — 



1 



1 



^Ci„(/) ^ 1.12 X 10-^Cin(/) 



2 I, Id 



(irA/)i/4 72 



A/ j 



1/4 

^ ) ■ (92) 



lyr 
T 



Of course, the reduction factor in the noise level is larger if the integration time is larger, 
and if we increase the bandwidth A/ over which a useful coincidence is possible. Note 
that /iQ^igw is quadratic in hc{f), so that an improvement in sensitivity by two orders of 
magnitudes in he means four orders of magnitude in /igfigw. 



5 Detectors in operation, under construction, or planned 

In this section we will discuss the characteristic of the detectors existing, under construction, 
or planned, and we will examine the sensitivities to stochastic backgrounds that can be 
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obtained using them as single detectors. The results will allow us to fully appreciate the 
importance of correlating two or more detectors, at least as far as stochastic backgrounds 
are concerned. 



5.1 The first generation of large interferometers: LIGO, VIRGO, GEO600, 
TAMA300, AIGO 

A great effort is being presently devoted to the construction of large scale interferometers. 
The Laser Interferometer Gravitational- Wave Observatory (LIGO) is being developed by 
an MIT-Caltech collaboration. LIGO consists of two widely separated interferometers, 
in Hanford, Washington, and in Livingston, Louisiana. The commissioning of the detectors 



will begin in 2000, and the first data run is expected to begin in 2002 |31|. The arms of 
both detectors will be 4 km long; in Hanford there will be also a second 2 km interferometer 
implemented in the same vacuum system. The sensitivity of the 4 km interferometers is 
shown in fig. |l|. The sensitivity is here given in meters per root Hz; to obtain hf, measured 
in I/VHz, one must divide by the arm length in meters, i.e. by 4000. The sensitivity in 
1/\/Hz is shown in fig. ^, next subsection, together with the sensitivity of advanced LIGO. 
The LIGO interferometer already has a 40m prototype at Caltech, which has been used to 
study sensitivity, optics, control, and even to do some work on data analysis 0]. 



A comparable interferometer is VIRGO |69, 256 , 53] which is being built by an Italian- 



French collaboration supported by INFN and CNRS, with 3 km arm length, and is presently 
under construction in Cascina, near Pisa. These experiments are carried out by very large 
collaborations, comparable to the collaborations of particle physics experiments. For in- 
stance, VIRGO involves about 200 people. Fig. |2| shows, in terms of hf, the expected 
sensitivity of the VIRGO interferometer. By the year 2000 it should be possible to have the 
first data from a 7 meters prototype, which will be used to test the suspensions, vacuum 
tube, etc. 

Interferometers are wide-band detectors, that will cover the region between a few Hertz 
up to approximately a few kHz. Figs. (|l|,§) show separately the main noise sources. At 
very low frequencies, / < 2 Hz for VIRGO and / < 40 Hz for LIGO, the seismic noise 
dominates and sets the lower limit to the frequency band. Above a few Hz, the VIRGO 
superattenuator reduces the seismic noise to a neglegible level. Then, up to a few hundreds 
Hz, thermal noise dominates, and finally the laser shot noise takes over. The various spikes 
are mechanical resonances due to the thermal noise: first, in fig. ^, the high frequency tail of 
the pendulum mode (in this region LIGO is still dominated by seismic noise), then a series 
of narrow resonances due to the violin modes of the wires; finally, the two rightmost spikes 
in the figures are the low frequency tail of the resonances due to internal modes of the two 
mirrors. All resonances appear in pairs, at frequencies close but non-degenerate, because 
the two mirrors have different masses. The width of these spikes is of the order of fractions 
of Hz. Note that the clustering of the various resonances in the kHz region is due to the 
logarithmic scale in figs, f^^- The resonances are actually evenly spaced and narrow, so 
that in this frequency range the relevant curve is mostly the one given by the shot noise. 



Under construction is also GEO600 [185], a collaboration between the Max-Planck- 



Institut fiir Quantenoptik in Garching and the University of Glasgow. It is being built 
near Hannover, with arms 600 meters in length. GEO600 is somewhat smaller than LIGO 



20 



Initial LIGO Noise Curves 
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Figure 1: The planned LIGO sensitivity curve, in meters/ v Hz (courtesy of the LIGO 
cohaboration) . 
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Frequency [Hz] 

Figure 2: The planned VIRGO sensitivity curve (from ref. [16]). 

and VIRGO, but will use techniques that will be important for the advanced detectors. In 
particular, the signal recycling of GEO600 provides an opportunity to change the spectral 
characteristics of the detector response, especially those due to the shot noise limitation, 
therefore in the high frequency range of the available bandwidth. By choosing low or 
high mirror reflectivities for the signal-recycling mirror, one can use the recycling either to 
distribute the improvement in the sensitivity on a wideband, or to improve it even further on 
a narrow band. This is the so-called narrow-banding of an interferometer. The sensitivity 
of GEO600 broad-band is shown in fig. ^, and with narrow-banding at 600 Hz in fig. ^. 

The frequency of maximum sensitivity is tunable to the desired value by shifting the 
signal-recycling mirror, and thus changing the resonance frequency of the signal-recycling 
cavity. In its tunable narrowband mode, GEO600 might well be able to set the most 
stringent limits of all for a while. 

TAMA300 (Japan) [163| is the other large interferometer. It is a five year project 
(1995-2000) that aims to develop the advanced techniques that will be needed for second 
generation experiments, and catch GWs that may occur by chance within out local group 
of galaxies. Its has 300 meters arm length. It is hoped that the project will evolve into the 
proposed Laser Gravitational Radiation Telescope (LGRT), which should be located near 
the SuperKamiokande detector. 

Furthermore, Australian scientists have joined forces to form the Australian Consortium 
for Interferometric Gravitational Astronomy (ACIGA) [ |227[ |. A design study and research 
with particular emphasis on an Australian detector is nearing completion. The detector, 
AIGO, will be located north of Perth. Its position in the southern emisphere will greatly 
increase the baseline of the worldwide array of detectors, and it will be close to the resonant 
bar NIOBE so that correlations can be performed. The detector will use sapphire optics 
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Figure 4: The planned GEO sensitivity curve, with narrow-banding at 600 Hz (courtesy of 
the GEO collaboration). 
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and sapphire test masses, a material that LIGO plans to use only in its advanced stage, see 
sect. IJ 

All these detectors will be used in a worldwide network to increase the sensitivity and 
the reliability of a detection. 

Using the values of /i/ given in the figures and the value 1/ F = 2.5 for interferometers, 
we see from eq. (54) that VIRGO, GEO600 or any of the two LIGOs , used as single 



detectors, can reach a minimum detectable value for h^Q.^ of order 10^^ or at most a few 
times 10~^, at / = lOOHz. Unfortunately, at this level current theoretical expectations 
exclude the possibility of a cosmological signal, as we will discuss in sect. and ^. As we 
will see, an interesting sensitivity level for H'qQ.^ should be at least of order 10~^. To reach 
such a level with a single interferometer we need, e.g., 

hf{f = 100 Hz) < 10^24 Hz-i/2 ^ (93) 

or 

hf{f = 1 kHz) < 3 X 10^2^ Hz"^/2 ^ (94) 

We see from the figures that such small values of /ij are very far from the sensitivity of first 
generation interferometers, and are in fact even well below the limitation due to quantum 
noise. Of course, one should stress that theoretical prejudices, however well founded, are no 
substitute for a real measurement, and that even a negative result at the level h^Q.^ ~ 10~^ 
would be interesting. 

5.2 Advanced LIGO 

It is important to have in mind that the interferometers discussed in the previous section 
are the first generation of large scale interferometers, and in this sense they represent really 
a pioneering effort. At the level of sensitivity that they will reach, they have no guaranteed 
source of detection. However, they will open the way to second generation interferometers, 
with much better sensitivity. GEO600 and TAMA300 are important for testing the tech- 
nique that will be needed for second generation experiments, while the larger interferometers 
VIRGO and LIGOs should evolve into second generation experiments. 

The LIGO collaboration has presented a recommended program for research and de- 



velopement, that will lead to the Advanced LIGO project |177]. The results are shown in 
fig. H, together with the changes in a number of parameters that are responsible for the 
various improvements shown in the figure. 

The improvement program is divided into three stages: LIGOII near term, LIGOII 
medium term and LIGOIII. 

Curve 1 in fig. |5| shows the sensitivity of LIGOI, while LIGOII near term is shown in 
curve 2. The changes leading to LIGOII near term should be incorporated by the end of 
2004, and are based on engeneering developents of existing technology or modest stretches 
from present day systems. They would result in significant improvements of the sensitivity, 
and are also necessary to gain full advantage of subsequent changes. In particular, a signif- 
icant reduction in thermal noise will be obtained using fused silica fibers in the test mass 
suspension; a moderate improvement in seismic isolation will move the seismic noise below 
this new thermal noise floor, and an increase in the laser power to ~ 100 watts results in a 
decrease of shot noise. 
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Figure 5: The planned advanced LIGO sensitivity curves, and the parameters relative to 
the various improvements, from ref. |177[ | (courtesy of the LIGO collaboration). 
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The goal of the medium term program are more ambitious technically, and it is estimated 
that they could be incorporated in LIGO by the end of 2006. They include signal recycling. 
As discussed above, this could be used either to obtain a broad-band redistribution of 
sensitivity (curve 3), or to improve the sensitivity in some narrow band (curve 4). 

Another possible improvement is the change of the test masses from fused silica to sap- 
phire or other similar materials. Sapphire has higher density and higher quality factor, and 
this should result in a significant reduction in thermal noise. For an optical configuration 
without recycling, the sensitivity is given by curve 5. With signal recycling, possible re- 
sponses are shown in curves 6 and 7. These changes require significant advances in material 
technology, and so are considered possible but ambitious for LIGOII. 

Fig. P does not show performances for LIGOIII detectors, which are expected to provide 
possibly another factor of 10 of improvement inhf, but of course are quite difficult to predict 
reliably at this stage. 

The overall improvement of LIGOII can be seen to be, depending on the frequency, one 
or two orders of magnitude in hf. This is quite impressive, since two order of magnitudes 
in /ij means four order of magnitudes in hQQg„{f) and therefore an extremely interesting 
sensitivity, even for a single detector, without correlations. 



5.3 The space interferometer LISA 



The space interferometer LISA 150 1 was proposed to the European Space Agency 

(ESA) in 1993, in the framework of ESA's long-term space science program Horizon 2000. 
The original proposal involved using laser interferometry beteween test masses in four drag- 
free spacecrafts placed in a Heliocentric orbit. In turn, this led to a proposal for a six space- 
crafts mission, which has been selected by the European Space Agency has a cornerstone 
mission in its future science program Horizon 2000 Plus. This implies that in principle the 
mission is approved and that funding for industrial studies and technology developement 
is provided right away. The launch year however depends on the availability of fundings. 
With reasonable estimates it is then expected that it will not be launched before 2017, and 
possibly as late as 2023. In Feb. 1997 the LISA team and ESA's Fundamental Physics 
Advisory Group proposed to carry out LISA in collaboration with NASA. If approved, this 
could make possible to launch it between 2005-2010. The design has also been somewhat 
simplified, with three drag-free spacecrafts. The spacecrafts will be in a Heliocentric orbit, 
at a distance of 1 AU from the Sun, 20 degrees behind the Earth. This design has lower 



costs, and it is likely to be adopted for future studies relevant to the project | 150 |. The 
mission is planned for two years, but it could last up to 10 years without exausting on-board 
supplies. 

LISA has three arms, first of all for redundancy. Thus, it can be thought of as two 
interferometers sharing a common arm. Of course, this means that the two interferometers 
will have a common noise. However, most signals are expected to have a signal-to-noise ratio 
so high that the noise will be neglegible. Then, the output from the two interferometers 
can be used to obtain extra informations on the polarization and direction of a GW. For 
the stochastic background, the third arm will help to discriminate backgrounds as those 



produced by binaries or by cosmological effects from anomalous instrumental noise |36] 



Goint into space, one is not limited anymore by seismic and gravity-gradient noises; 
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LISA could then explore the very low frequency domain, 10~^ Hz < / < 1 Hz. At the 
same time, there is also the possibility of a very long path length (the mirrors will be 
freely floating into the spacecrafts at distances of 5 x 10^ km from each other!), so that 
the requirements on the position measurment noise can be relaxed. The goal is to reach a 



strain sensitivity |36| 



hf = Ax 10"^^ Hz"^/2 (95) 

at / = 1 mHz. At this level, one expects first of all signals from galactic binary sources, 
extra-galactic supermassive black holes binaries and super-massive black hole formation. 

Concerning the stochastic background, eq. (^) shows clearly the advantage of going 
to low frequency: the factor in eq. (|5^) gives an extremely low value for the minimum 
detectable value of /igOgw. The sensitivity of eq. (^) would correspond to 

/io^gw(/ = 1 mHz) ~ 1 X 10^^^ . (96) 

Since LISA cannot be correlated with any other detector, to have some confidence in the re- 
sult it is necessary to have a SNR sufficiently large. Certainly one cannot work at SNR=1.65, 
as we will do when considering the correlation between two detectors. The standard choice 
made by the LISA collaboration is SNR=5, and eq. (|9^ ) refers to this choice. Note that 
the minimum detectable value of h^Qgw is proportional to (SNR)^, since this SNR refers to 
the amplitude, and /igJ7gw(/) ~ ^c(/)- Furthermore, eq. (|96|) takes into account the angles 
between the arms, a = 60°, and the effect of the motion of LISA, which together results 
in a loss of sensitivity by a factor approximately equal to 1/V5. Eq. (|9^ ) shows that LISA 
could reach a truly remarkable sensitivity in /igilgw. 

LISA will have its best sensitivity between 3 and 30 mHz. Above 30 mHz, the sensitivity 
degrades because the wavelength of the GW becomes shorter than twice the arm-length of 
5 X 10^ km. At low frequencies, instead, the noise curve rises because of spurious forces on 
the test masses. At some frequency below 0.1 mHz the accelerometer noise will increase 
rapidly, and the instrumental uncertainty would increase even more rapidly with decreasing 
frequency, setting a lower limit to the frequency band of LISA. 

Below a few mHz it is expected also a stochastic background due to compact white- 
dwarf binaries, that could cover a cosmological background. The sensitivity curve of LISA 
to a stochastic background, together with the estimated white-dwarf binaries background, 
is shown in fig. On the vertical axis is shown /ir.m.s. = hc{f,Af = 3 x 10~^Hz), defined 



in eq. (|2^, which is written in |36] in the equivalent form 



/.,(/,A/ = 3xlO-SHz) = 5.5x10-22 fi^V^Vi^V^V ^ -) . (97) 

' VlO-V \ f J V75kms-iMpc-V ^ ^ 

Two lines of constant Ogw, equal to 10^^^ and 10^*^, are also shown in the figure. 

To compare with interferometers, fig. |^ shows the sensitivity to the GW amplitude 
for both LISA and the advanced LIGO, together with a number of signals expected from 
astrophysical sources. It is apparent from the figure that ground-based and space-borne 
interferometers are complementary, and together can cover a large range of frequencies, 
and, in case of a detection of a cosmological signal, together they can give crucial spectral 
informations. 

The sensitivities shown are presented by the collaboration as conservatives because 
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Figure 6: The sensitivity of LISA to a stochastic background of GWs after one year of 
observation, and SNR=5 (from ref. |36]; courtesy of the LISA cohaboration). 
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Figure 7: The sensitivity of LISA for periodic sources, and of advanced LIGO for burst 
sources, together with some expected astrophysical signals (from ref. [^6|; courtesy of the 
LISA collaboration). 
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1. The errors have been calculated realistically, including all substantial error sources 
that have been thought of since early studies of drag-free systems, and since the first 
was flown over 25 years ago, and in most cases (except shot noise) the error allowance 
is considerably larger than the expected size of the error and is more likely an upper 
bound. 

2. LISA is likely to have a significantly longer lifetime than one year, which is the value 
used to compute these sensitivities. 

3. The sensitivity shown refers only to one interferometer. Using three arms could in- 
crease the SNR by perhaps 20% . 

An interesting feature of LISA is that, as it rotates in its orbit, its sensitivity to different 
directions changes. Then, LISA can test the isotropy of a stochastic background. This can 
be quite important to separate a cosmological signal from a stochastic background of galactic 
origin, which is likely to be concentrated on the galactic plane. By comparing two 3-month 
stretches of data, LISA should have no difficulty in identifying this effect. Furthermore, 
if the mission lasts lOyr, LISA will be close to the sensitivity needed to detect the dipole 
anisotropy in a cosmological background due to the motion of the solar system. If the GW 
background turned out not to have the same dipole anisotropy as the cosmic microwave 
background, one would have found evidence for anisotropic cosmological models. 



5.4 Resonant bars: NAUTILUS, EXPLORER, AURIGA, ALLEGRO, 
NIOBE 



Cryogenic resonant antennas have been taking date since 1990 (see e.g. ref. |215] for review 



Fig. ^ shows the sensitivity of NAUTILUS [17, [l^, |7^, the resonant bar located in 



Frascati, near Rome. It is an ultracryogenic detector, operating at a temperature T = 0.1 
K. This figure is based on a two hour run, but the behaviour of the apparatus is by now 
quite stationary over a few days, with a duty cycle limited to 85% by cryogenic operations. 
Compared to the interferometers, we see that bars are narrow-band detectors, and work at 
two resonances. The bandwidth is limited basically by the noise in the amplifier. The value 
of the resonance can be slightly tuned, at the level of a few Hz, working on the electronics. 
In the figure, the resonances are approximately at 907 Hz and 922 Hz, with half-height 
bandwiths of about 1 Hz. At these frequencies the strain sensitivity is 

/i/ ~ 5 X 10^22 , (98) 

about a factor of 5 higher than the sensitivity of the interferometers at the same frequencies. 

With a thermodynamic temperature of 0.1 K, the experimental data are in very good 
agreement with a gaussian distribution, whose variance gives an effective temperature T^s = 
4.1 mK. The effective temperature T^s is related to the thermodynamic temperature T by 
Tgff ~ 4Tr^/^, where F is inversely proportional to the quality factor Q of the resonant 
mode, so that F <C 1. 

With improvements in the electronics, the collaboration plans to reach in a few years 
the target sensitivity 

/i/ ~ 8.6 X lO^^^Hz"^/^ (99) 
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Figure 8: The sensitivity curve of NAUTILUS (courtesy of the NAUTILUS collaboration) 



and TcfT WfiK, over a bandwidth of 10 Hz 

The cryogenic resonant bar EXPLORER P08 ] is in operation at Cern since 1990, cooled 
at a temperature T = 2.9 K, and has taken data at approximately the same frequencies as 
NAUTILUS, 905 and 921 Hz; its sensitivity is approximately the same as NAUTILUS. 

Similar sensitivities are obtained by ALLEGRO, (Louisiana) 



|193| , |14|. ALLEGRO has 
been taking data since 1991. The resonances are at 897 and 920 Hz. 

The resonant bar NIOBE ES, E9, 144, |145|1 is located in Perth, Australia, and it has 



operated since 1993. The resonant frequencies are at 694.6 Hz and 713 Hz. It is made of 
niobium, instead of aluminium as the other bars, since it has a higher mechanical quality 
factor. This has allowed to reach a noise temperature of less then 2mK, while the typical 
noise temperature is around 3mK. Improvement in the transducer system should allow to 
lower Teff to a few microkelvin, in a bandwidth of 70 Hz. This would allow to reach a 
sensitivity comparable to interferometers over a large bandwidth, approximately 650-750 
Hz p5|l . 

1911 is located in Legnaro, near Padua, Italy. It began full operations 

7 mK. The 



218 



AURIGA 

in 1997. The bar is cooled at 0.2K, and the effective temperature is Tcfj 
resonances are at 911 and 929 Hz. The spectral sensitivity reached |21£] is similar to that 
of NAUTILUS. 

All these bars are designed to operate in a well coordinated way, in order to improve 
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the chances of rehable detection^ 

From eq. ( |99|) we see that, without performing correlations between different detectors, 
the sensitivity that one could get from resonant bars is between at most hQQ'^[f) ~ 
0(1) and /iqO™"(/) ~ 0(10) and therefore quite far from a level where one can expect a 
cosmological signal. In sect. ^ we will discuss the level that can be reached with correlations 
between bars, and between a bar and an interferometer. 



5.5 Some projects at a preliminary stage 



A number of other projects for gravitational wave detection have been put forward and are 
presently at the stage of testing prototypes. 

An especially interesting project involves resonant mass detectors with truncated icosa- 

|]. A spherical detector, in fact, has a larger mass 



hedron geometry ^5^, |l9|, |9|, |g, |2 



compared to a bar with the same resonance frequency, and therefore a larger cross section. 
Furthermore, it has the informations about the directions and polarization of a GW that 
could only be obtained with 5 different resonant bars. The problem of attaching mechanical 
resonators to the detector suggests a truncated icosahedral geometry, rather than a sphere. 
A prototype called TIGA (Truncated Icosahedral Gravitational Wave Antenna) has 



been built at the Lousiana State University [198]. It has its first resonant mode at 3.2 



kHz. The Rome group has started the SFERA project [25|. A working group has been 
formed to to carry out studies and measurements in order to define a project of a large 
spherical detector, 40 to 100 tons of mass, competitive with large interferometers but with 
complementary features. A correlation between an interferometer and a sphere could in 
particular have quite interesting sensitivities for the stochastic background. 

These spheres could reach approximately one order of magnitude better in the sensitivity 
182 1 compared to resonant bars. Furthermore, they have different channels, see 
that are sensitive to different spin content |11^, p 



) and 



eqs 

allow to discriminate between different theories of gravity, separating for instance the effects 



157]. This would 



of scalar GWs, which could originate from a Brans-Dicke theory |255, 159, 182, 54, 42, 



206, 189]. Scalars fields originating from string theory, as the dilaton and the moduli of 
compactifications, are however a much more difficult target: to manifest themselves as 
coherent scalar GWs they should be extremely light, m < 10"^'^ eV J189U . Actually, there 
is the possibility to keep the dilaton and moduli fields extremely light or even massless, 
with a mechanism that has been proposed by Damour and Polyakov ]^^. In fact, assuming 
some form of universality in the string loop corrections, it is possible to stabilize a massless 
dilaton during the cosmological evolution, at a value where it is essentially decoupled from 
the matter sector. In this case, however, the dilaton becomes decoupled also from the 
detector, since the dimensionless coupling of the dilaton to matter (a in the notation of ]p4| ) 
is smaller than 10""^ (see also |9^). Such a dilaton would then be unobservable at VIRGO, 
although it could still produce a number of small deviations from General Relativity which 
might in principle be observable improving by several orders of magnitude the experimental 



tests of the equivalence principle ]9C]. 



^There is an International Gravitational Event Collaboration, an agreement between the bar detector 
experiments presently in operation, signed at CERN on July 4, 1997. 
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Recently, hollow spheres have been proposed in ref . ||78| . The theoretical study suggests 
a very interesting value for the strain sensitivity, even of order 



hf ~ afew X 10~^^ . (100) 

The resonance frequency, depending on the material used and other parameters, can be 
between approximately 200 Hz and 1-2 kHz, and the bandwidth can be of order 20 Hz. 



Another idea which is discussed is an array of resonant masses |lig| , 120, 111|, each 
with a different frequency /, and a bandwidth Af //lO, which together would cover the 
region from below the kHz up to a few kHz. 

Although not of immediate applicability to GW search, we finally mention the proposal 
of using two coupled superconducting microwave cavities to detect very small displacements. 



The idea goes back to the works |211, ^12| , |213| , |153| . These detectors have actually been 



constructed in ref. [p21[| , using two coupled cavities with frequencies of order of lOGHz; the 
coupling induces a shift in the levels of order of 1 MHz, and transitions between these levels 
could in principle be induced by GWs. However, eq. (|5^ ) shows clearly the problem with 
working at such a high frequency, / ~ 1 MHz: to reach an interesting value for /iQilgw(/) 
of order 10~ , one needs hf ~ 10~^° Hz'^/^ very far from present technology. There is 
however the proposal |40| to repeat the experiment with improved sensitivity, so to reach 



hf ^6x IQ-^^ Hz~^/2 at / = 1 MHz, and, if this should work, then one would try to lower 
the resonance frequency down to, possibly, a few kHz. This second step would mean to 
move toward much larger cavities, say of the kind used at LEP. 



6 Sensitivity of various two-detectors correlations 

In the previous section we have seen that with a single detector (and with the exception of 
LISA) we cannot reach sensitivities interesting for cosmological backgrounds of GWs. On 
the other hand, we have found in sects. (^^) and ([4.4| ) that many orders of magnitudes can 
be gained correlating two detectors. In this section we examine the sensitivities for various 
two-detector correlations. 

It is also important to stress that performing multiple detector correlations is crucial 
not only for improving the sensitivity, but also because in a GW detector there are in 
general noises which are non-gaussian, and can only be eliminated correlating two or more 
detectors. An example of a non-gaussian noise in an interferometer is the creep, i.e. the 
sudden energy release in the superattenuator chain, and most importantly in the wires 
holding the mirrors, and due to accumulated internal stresses in the material. Indeed, 
it is feared that this and similar non-stochastic noises might turn out to be the effective 



sensitivity limit of interferometers, possibly exceeding the design sensitivity \98]. Some 
sources of non-stochastic noises have been identified, and there are technique for neutralizing 
or minimizing them. However, it is impossible to guarantee that all possible sources of non- 
stochastic noises have properly been taken into account, and only the cross correlation 
between different detectors can provide a reliable GW detection. 
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6.1 An ideal two-interferometers correlation 



We now discuss the sensitivities that could be obtained correlating the major interferome- 
ters.| 

In order to understand what is the best result that could be obtained with present inter- 
ferometer technology, we first consider the sensitivity that could be obtained if one of the 
interferometers under constructions were correlated with a second identical interferometer 
located at a few tens of kilometers from the first, and with the same orientation. This dis- 
tance would be optimal from the point of view of the stochastic background, since it should 
be sufficient to decorrelate local noises like, e.g., the seismic noise and local electromagnetic 
disturbances, but still the two interferometers would be close enough so that the overlap 
reduction function does not cut off the high frequency range. For this exercise we use the 
data for the sensitivity of VIRGO, fig. |2|. 

Let us first give a rough estimate of the sensitivity using hc{f), hn{f). Prom fig. || we 
see that we can take, for our estimate, ^/ ~ 10~^^Hz~^/'^ over a bandwidth A/ ~ 1 kHz. 
Using T = 1 yr, eq. ( |88[ ) gives 

.-24 f f h ' 



KU) ~ 4.5 X 10-^4 '—-^ . (101) 

^ ' VlOOHzy \^10-22Hz-i/V ^ ' 

We require for instance SNR=1.65 (this corresponds to 90% confidence level; a more precise 
discussion of the statistical significance, including the effect of the false alarm rate can be 
found in ref. [ff^]). Then the minimum detectable value of hc{f) is h^^^{f) = 1.65/i„(/), 
and from eq. (pj|) we get an estimate for the minimum detectable value of /ioi^gw(/), 



K^fwif) ~ 3 X lO-M — ^ ^—7^ . (102) 

° ^^^■'^ VlOOHzy yiO-22Hz~i/V ^ ^ 

This suggests that correlating two VIRGO interferometers we can detect a relic spectrum 
with /igrig„(100Hz) ~ 3 X 10""^ at SNR=1.65, or 1 x 10""^ at SNR=1. Compared to the 
case of a single interferometer with SNR=1, eq. (|54[), we gain five orders of magnitude. As 
already discussed, to obtain a precise numerical value one must however resort to eq. ([79|). 
This involves an integral over all frequencies, (that replaces the somewhat arbitrary choice of 
A/ made above) and depends on the functional form of hQQg„{f). If for instance /iQJ7gw(/) 
is independent of the frequency, using the numerical values of hf plotted in fig. ^ and 
performing the numerical integral, we have found for the minimum detectable value of 
/)2o 

.2omin^o.,n-7^SNR^Vlyr^^/^ 



h',n'^-^2xW-^ ^-^^ ^-jr) (/i^f^g„(/) = const.) . (103) 

This number is quite consistent with the approximate estimate (|102|), and with the value 



2 X 10 ' reported in ref. Stretching the parameters to SNR=1 (68% c.l.) and T = 4 



® Correlations between two interferometers have already been carried out using prototypes operated 
by the groups in Glasgow and at the Max Planck Institute for Quantum Optics, with an effective coicident 



observing period of 62 hours [205]. Although the sensitivity of course is not yet significant, they demonstrate 



the possibility of making long-term coincident observations with interferometers. 
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years, the value goes down at (3 — 4) x 10^®. This might be considered an absolute (and 
quite optimistic) upper bound on the capabilities of first-generation experiments. 

It is interesting to note that the main contribution to the integral comes from the region 
/ < 100 Hz. In fact, neglecting the contribution to the integral of the region / > 100 Hz, the 
result for /ipfi™" changes only by approximately 2%. Also, the lower part of the accessible 
frequency range is not crucial. Restricting for instance to the region 20 Hz < f < 200 Hz, 
the sensitivity on /igOgw degrades by less than 1%, while restricting to the region 30 Hz 
< / < 100 Hz, the sensitivity on /ig^^gw degrades by approximately 10%. Then, from fig. ^ 
we conclude that by far the most important source of noise for the measurement of a flat 
stochastic background is the thermal noise. In particular, the sensitivity to a flat stochastic 
background is limited basically by the mirror thermal noise, which dominates in the region 
40 Hz</<200 Hz, while the pendulum thermal noise dominates below approximately 40 
Hz. 

The sensitivity depends however on the functional form of 0,g^{f). Suppose for instance 
that in the VIRGO frequency band we can approximate the signal as 

J^gw(/) = l^gw(lkHz) (^j^y . (104) 

For a = 1 we find that the spectrum is detectable at SNR=1.65 if /lo^^gw(lkHz) ~ 3.6x 10~^. 
For a = —1 we find (taking / = 5Hz as lower limit in the integration) /igOgw(lkHz) ~ 
6 X 10^^. Note however that in this case, since q < 0, the spectrum is peaked at low 
frequencies, and /iQ$7g„(5Hz) ~ 1 x 10^^. So, both for increasing or decreasing spectra, to 
be detectable /ig^gw must have a peak value, within the VIRGO band, of order a few x 10~^ 
in the case a = ±1, while a constant spectrum can be detected at the level 2 x 10"'''. Clearly, 
for detecting increasing (decreasing) spectra, the upper (lower) part of the frequency band 
becomes more important, and this is the reason why the sensitivity degrades compared to 
fiat spectra, since for increasing or decreasing spectra the maximum of the signal is at the 
edges of the accessible frequency band, where the interferometer sensitivity is worse. 



6.2 LIGO-LIGO 

Let us now see what can be done with existing interferometers. The two LIGO detectors are 
under construction at a large distance from each other, d ~ 3000 km. This choice optimizes 
the possibility of detecting the direction of arrival of GWs from astrophysical sources, but 
it is not optimal from the point of view of the stochastic background, since the overlap 
reduction function cuts off the integrand in eq. ( [7^ ) at a frequency of the order of l/{2TTd). 
The overlap function 7(/) for the LIGO-LIGO correlation has been computed in ref. [|113| . 
and it has its first zero at / ~ 64 Hz. Furthermore, the arms of the two detectors are not 
exactly parallel, and therefore |7(0)| = 0.89 rather than 1. 

The sensitivity to a stochastic background for the LIGO-LIGO correlation has been 



computed in refs. |200, 74, 112, E, O]. The result, as we have discussed, depends on the 



functional form of /iQ$7gw. For Qg„(f) independent of /, the minimum detectable value is 

/lo^gw =^ 5 X 10"^ (105) 

for the initial LIGO. However, we have seen that second generation interferometers could 
result in much better sensitivities. The sensitivity of the correlation between two advanced 
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correlation 


gw 


LIGO-WA*LIGO-LA 


5 X 10"*^ 


VIRGO*LIGO-LA 


4 X 10-*^ 


VIRGO*LIGO-WA 


5 X 10-*^ 


VIRGO*GEO600 


5.6 X 10"^ 


VIRGO*TAMA300 


1 X 10"^ 



Table 1: The sensitivity of various two- interferometers correlation. 



LIGO is estimated in ref. |^] to be 

/igJlgw ^ 5 X 10"" , (106) 

which is an extremely interesting level. These numbers are given at 90% c.l. in and a 
detailed analysis of the statistical significance is given in |12|. 

6.3 VIRGO-LIGO, VIRGO-GEO, VIRGO-TAMA 

In Fig. P we show the overlap reduction functions for the correlation of VIRGO with the 
other major interferometers. Using these functions, one can compute numerically the inte- 
gral in eq. ([79| ) and obtain the minimum detectable value of /iQ$7g„. These values are shown 
in Table 1 (from ref. [^]), together with the value for LIGO-LIGO computed in ref. All 
these numbers are at 90% confidence level. 

We see that the correlation between VIRGO and any of the two LIGO is suppressed 
by the overlap reduction function above, say, 30-40 Hz. In the case of the VIRGO-GEO 
correlation, instead, 7(/) cuts the integrand only above say 200 Hz. However, the sensitivity 
of GEO is below 200 Hz is lower then LIGO, so that the result are basically the same as 
for VIRGO-LIGO, see Table 1. 

To obtain the precise numbers for the sensitivity, of course one has to perform the 
integral in eq. (^). However, it is easy to have an understanding of the numbers that 
come out. As we have found in sect. 6.1, the minimum detectable value of h^Qg^, using 
two identical VIRGO detectors, degrades only by less than 10% if, from the full VIRGO 
bandwidth, we restrict to the region 20 Hz < / < 100 Hz; in this region |7(/)| for the 



VIRGO-GEO correlation is constant to a good accuracy, and of order 0.1. From eqs. ( |79[) 
and (|l^) we see that the minimum detectable value of /igf^gw scales with 7(/) as |7(/)|~^. 
Therefore for the VIRGO-GEO correlation we get about a factor of 10 worse than the ideal 
result ( [1031) . Comparing figs. |2| and ^, and recalling that /igf^™"^ ~ /i^, we see that we lose 
about another factor of 3 in /igfigw compared to an ideal VIRGO- VIRGO correlation, so 
that one gets the estimate 

hln^:- ~ 6 X 10-6 ' ^'^ {hln,M = const.) , (107) 



quite in agreement with the result of the more accurate computation shown in Table 1. 

The situation with the VIRGO-TAMA300 correlation is instead worse, as can be seen 
from the overlap reduction function and from the value in Table 1. 
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correlation 




VIRGO*AURIGA 


4 X 10-"^ 


VIRGO*NAUTILUS 


7 X lO-'* 


AURIGA*NAUTILUS 


5 X 10""^ 



Table 2: The sensitivity of the correlation of VIRGO with resonant bars, and of two resonant 
bars (from ref. [^9|), but setting the confidence level to 90% ). 



6.4 VIRGO-Resonant mass and Resonant mass-Resonant mass 

The correlation between two resonant bars and between a bar and an interferometer has 
been considered in refs. 

Fig. ^ shows the overlap reduction functions for the correlation between VIRGO and 
one of the three resonant bars NAUTILUS, EXPLORER, AURIGA. These overlap reduction 
functions have been computed assuming that the bars have been reoriented so to achieved 
the maximum correlation with VIRGO, which is technically feasible. 

One should also note that there is the danger that one of the spikes due to the violin 
modes in the VIRGO sensitivity curve comes close to a resonant frequency of the bar; 
surprisingly, this apparent unlikely event is just what happens with the data used to draw 
fig. in fact, according to these data there is a pair of violin modes at 922.6 Hz and 977.2 Hz. 
The first one happens to fall exactly at the resonance at 922 Hz of NAUTILUS! However, 
these data for the violin modes are not yet final. For instance, the VIRGO collaboration 
is presently considering the possibility of using silica instead of steel for the wires, which 
would change the position of the resonances. Furthermore, the resonance frequency of the 
bars can be tuned within a few Hz with the electronics; this would be quite sufficient, since 
the violine modes have a very high Q and so are much narrower than one Hz. 

The minumum detectable values for /igf^gw for some bar-bar and bar-interferometer 
correlations are given in Table 2 (from ref. |^9[), for one year of observation and 90% 
confidence level (SNR=1.65). 

A three detectors correlation AURIGA-NAUTILUS-VIRGO, with present orientations, 
would reach 

/iQ^gw ^ 3 X 10""^ , (108) 
at 90% C.I., while with optimal orientation, 

/igjlgw ~ 1.6 X 10~^. (109) 

Although the improvement in sensitivity in a bar-bar-interferometer correlation is not large 
compared to a bar-bar or bar-interferometer correlation, a three detectors correlation would 
be important in ruling out spurious effects |257|] . 

In the case of correlations involving resonant bars, obviously there is no issue of optimal 
filtering, since they are narrow-band detector, and the sensitivity does not depend on the 
shape of the spectrum. As we have seen in sect. |5|, using resonant optical techniques, it is 
possible to improve the sensitivity of interferometers at special values of the frequency, at 



257, 83, 21, H, 23 
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Figure 10: The overlap reduction function 7(/) for the correlation of VIRGO with NAU- 
TILUS (sohd hue), AURIGA (dashed line) and EXPLORER (dotted hue) in the case in 
which the bars have been reoriented so to be aligned with VIRGO (from ref. [p9|]). 
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the expense of their broad-band sensitivity. Since bars have a narrow band anyway, narrow- 
banding the interferometer improves the sensitivity of a bar-interferometer correlation by 
about one order of magnitude [^]. Thus, the Umit of bar-bar-interferometer correlation, 
with narrow banding of the interferometer, is of order /igf^gw ~ a few x 10~^. 

Cross-correlation experiments have already been performed using NAUTILUS and EX- 
PLORER 1^^. The bar are oriented so that they are parallel, and then the overlap reduction 
function r(/) results in a reduction of sensitivity of about a factor of 6, compared to ideal 
same site detectors. The detectors are tuned at the same resonance frequency / = 907.2 
Hz, with an overlapping band ±0.05 Hz, and the overlapped data cover a period of approx- 
imately 12 hours. Of course, with such a short coincidence time, the bound on /igl^gw is not 
yet significant. Running for one year, it is expected to reach /igfigw < 1- Cross-correlations, 
searching for bursts, have also been done between ALLEGRO and EXPLORER 

It is also important to observe that, at the level of sensitivity that we are discussing, the 



effect of cosmic rays on the resonant detectors can be relevant [13, 82, 194, 26, 216|. For a 



stochastic background, cosmic rays at sea level degrade the power spectrum sensitivity of 



a single resonant mass detector by a factor (20mK + Teff)/Teff |216| . For the coincidence 
between two resonant masses it would then be useful to place one of the two detectors 
underground (not much is gained placing both detectors underground). In the coincidence 
between a resonant mass and an interferometer, instead, this is not necessary, since the 
interferometer is much less sensitive to cosmic rays. 

While resonant bars have been taking data for years, spherical detectors are at the 
moment still at the stage of theoretical studies (although prototypes might be built in the 
near future), but could reach extremely interesting sensitivities. In particular, two spheres 
with a diameter of 3 meters, made of A15056, and located at the same site, could reach a 



sensitivity /ig^gw ~ 4 x 10 [257|. This figure improves using a more dense material or 



increasing the sphere diameter, but it might be difficult to build a heavier sphere. Another 



very promising possibility is given by hollow spheres [78|. The theoretical studies of ref. |]7 



suggest that correlating two hollow spheres one could reach the 

>'l2ok) 10-^^4-/0 (^r) —1 ■ <™> 



With the value of /i/ and A/ suggested by ref. [^], it could be possible to reach an extremely 



interesting value, /igi^gw ~ 10 ^ 
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7 Bounds on hlQg^{f) 

In this section we discuss various experimental bounds on /igr2gw(/). We will be interested 
not only on the bounds at values of / where interferometers or resonant bars can operate, 
but also at all possible frequencies. The reason will become apparent when we will discuss 
the spectra from various specific cosmological mechanisms for the production of relic GWs. 
These spectra depend of course on the parameters of the cosmological model. Often the 
frequency dependence is, to a first approximation, completely determined, but the overall 
value of /igOgw depends on some parameters of the model. In some case, and especially for 
the amplification of vacuum fluctuations (sect. |9.1D these spectra extend over a huge range 
of frequencies, ranging from frequencies as small as 10~^^ Hz (corresponding to wavelength 
of the order of the present Hubble radius of the Universe) up to possibly the GHz region. 
It is therefore important to see what are the experimental constraint, at any frequency, on 
/igOgw, since they automatically imply bounds on the parameters of the model that enter 
in the spectrum, and therefore on its the value at frequencies accessible to interferometers 
or resonant masses. 



The various limits discussed in this section are summarized in fig. 11 
7.1 The nucleosynthesis bound 

Nucleosynthesis successfully predicts the primordial abundances of deuterium, ^He, '^He 
and ^Li in terms of one cosmological parameter r], the baryon to photon ratio. In the 
prediction enter also parameters of the underlying particle theory, which are therefore con- 
strained in order not to spoil the agreement. In particular, the prediction is sensitive to 
the effective number of species at time of nucleosynthesis, = g{T ~ MeV). With some 
simplifications, the dependence on g^, can be understood as follows. A crucial parameter 
in the computations of nucleosynthesis is the ratio of the number density of neutrons, n^, 
to the number density of protons, n^. As long as thermal equilibrium is mantained we 
have (for non-relativistic nucleons, as appropriate at T ~ MeV, when nucleosynthesis takes 
place) rin/up = exp{—Q/T) where Q = rrin — mrip ~ 1.3 MeV. Equilibrium is mantained by 
the process pc <^ Tiiy y with, width. Tpe >niy-) 

as long as Tpe^n^ > H. When the rate drops 
below the Hubble constant H, the process cannot compete anymore with the expansion of 
the Universe and, apart from occasional weak processes, dominated by the decay of free 
neutrons, the ratio rin/np remains frozen at the value exp(— Q/Tf), where Tf is the value 
of the temperature at time of freeze-out. This number therefore determines the density 
of neutrons available for nucleosynthesis, and since practically all neutrons available will 
eventually form ^He, the final primordial abundance of ^He is exponentially sensitive to 
the freeze-out temperature Tf. Let us take for simplicity Tp^-^nu ~ CpT^ (which is really 
appropriate only in the limit T ^ Q). The Hubble constant is given by = (87r/3)G/9, 
where p includes all form of energy density at time of nucleosynthesis, and therefore also 
the contribution of primordial GWs. As usual, it is convenient to write the total energy 
density p in terms of 5*, see eq. ( |152 ), as p = {tt"^ / 30) g^T^. We recall that, for gravitons. 



the quantity Tj entering eq. ( |152| ) is defined by />gw = 2{7r'^ / 30)Tf , and this does not imply 
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Figure 11: The various bounds on hQilg^fif) discussed in the text. 
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a thermal spectrum. Then the freeze-out temperature Tf is determined by the condition 



1/2 rr.2 

90 / Mpi 



alTfJ'-^Y 4!- (HI) 



1 /6 

This shows that T{ gj , at least with the approximation that we used for Tp^^nu- A 
large energy density in relic gravitons gives a large contribution to the total density p and 
therefore to g^:. This results in a larger freeze-out temperature, more available neutrons and 



then in overproduction of ^He. This is the idea behind the nucleosynthesis bound [233|. 
More precisely, since the density of ^He increases also with the baryon to photon ratio rj, we 
could compensate an increase in g^, with a decrease in rj, and therefore we also need a lower 
limit on r], which is provided by the comparison with the abundance of deuterium and ^He. 

Rather than g^^, it is often used an 'effective number of neutrino species' Ni, defined 
as follows. In the standard model, at T ~ a few MeV, the active degrees of freedom are 
the photon, e^, neutrinos and antineutrinos, and they have the same temperature, Tj = T. 
Then, for Ni, families of light neutrinos, g>f{Ny) = 2 + (7/8) (4 + 2Ny)^ where the factor of 
2 comes from the two elicity states of the photon, 4 from in the two elicity states, and 
2iVi/ counts the N^, neutrinos and the Ny antineutrinos, each with their single elicity state; 
the factor (7/8) holds for fermions. In the Standard Model with N^, = 3, g^, = 43/4. So we 
can define an 'effective number of neutrino species' Ni, from 

t'+ E a^^y + l E m(^X = 2 + l{A + 2N.), (112) 

or 



4 . . u \T J 8 . V . V T 

2=extra bosons i=extraiermions 



E 9.(^y + l E 9^(^y = l{N.-3). (113) 

i=extra bosons i=extra fermions 

One extra species of light neutrino, at the same temperature as the photons, would con- 
tribute one unit to N^, but all species, weighted with their energy density, contribute to 
Ni,, which of course in general is not an integer. For i = gravitons, we have gi = 2 and 
{Ti/T)^ = pg^/p^, where p^ = 2(7r^/30)r^ is the photon energy density. If gravitational 
waves give the only extra contribution to N^, compared to the standard model with = 3, 
then 



=extra bosons 



and therefore 




..(|) =2^. (114) 



(115) 



NS 

where the subscript NS reminds that this equality holds at time of nucleosynthesis. If more 
extra species, not included in the standard model, contribute to g*{N^), then the equal 
sign in the above equation is replaced by lower or equal.[^ The same happens if there is a 

^To compare with eq. (56) of ref. note that while we use p-y, in ref. the bound is written in terms 
of the total energy density in radiation at time of nucleosynthesis, which includes also the contribution of 
e*, neutrinos and antineutrinos, (prad)NS = [1 + (7/8)(2 + 3)](p7)ns- 
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contribution from any other form of energy present at the time of nucleosynthesis and not 
included in the energy density of radiation, like, e.g., primordial black holes. 

To obtain a bound on the energy density at the present time, we note that from the time 
of nucleosynthesis to the present time /jgw scaled as while the photon temperature 

evolved following gs{T)T^a? = const, (i.e. constant entropy) and pj 1/ {a^gg^). 

Therefore 

(p^\ ^(P^\ ( 9s{To) y/^_/pgw\ / 3.913 

where the subscript zero denotes present time. Therefore we get the nucleosynthesis bound 
at the present time, 

j^j < 0.227(iV^ - 3) . (117) 
yp-y Jo 

Of course this bound holds only for GWs that were already produced at time of nucle- 
osynthesis (T ~ MeV, t ~ sec). It does not apply to any stochastic background pro- 
duced later, like backgrounds of astrophysical origin. Note that this is a bound on the 
total energy density in gravitational waves, integrated over all frequencies. Writing /?g„ = 
/ d{log f) dpgw/dlog/, multiplying both pg^ and p.y in eq. ( |117| ) by h^/pc and inserting the 
numerical value h^p^/ pc — 2.481 x 10^^ [l209|] , we get 

/•f=oo 

/ d(log/) /igf]gw(/) < 5.6 X 10-6(iV, - 3) . (118) 
Jf=o 

The bound on Ni, from nucleosynthesis is subject to various systematic errors in the analysis, 
which have to do mainly with the issues of how much of the observed '^He abundance is 
of primordial origin, and of the nuclear processing of '^He in stars, and as a consequence 
over the last five years have been quoted limits on Ny ranging from 3.04 to around 5. The 



situation has been reviewed in ref. [89|. The conclusions of ref. |89] is that, until current 
astrophysical uncertainties are clarified, N^, < 4 is a conservative limit. Using extreme 
assumptions, a meaningful limit Ni, < 5 still exists, showing the robustness of the argument. 
Correspondingly, the right-hand side of eq. ( |118D is, conservatively, of order 5 x 10~^ and 
anyway cannot exceed 10~^. 

If the integral cannot exceed these values, also its positive definite integrand hQQg^{f) 
cannot exceed it over an appreciable interval of frequencies, A log / ~ 1. One might still 
have, in principle, a very narrow peak in /iQilgw(/) at some frequency /, with a peak value 
larger than say 10~^, while still its contribution to the integral could be small enough. 
But, apart from the fact that such a behaviour seems rather implausible, or at least is 
not suggested by any cosmological mechanism, in the detection at broadband detectors like 
VIRGO, one must also take into account that, even if we gain in the height of the signal, 
this is partially compensated by the fact that we lose because of the reduction of the useful 



frequency band A/, see eq. (|78D . 



It is sometimes useful to think in terms of the dimensionless amplitude /ic(/) instead of 
/igilgw. From eq. (|T9|), we see that a bound hQQg^{f) < 5 x 10~^, valid independently of 
the frequency /, translates into a frequency dependent bound on hdf), 

Kif) < 2.82 X 10-21 (^i^) , (119) 
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so that, for instance, hdf) < 10 around 200-300 Hz, where ground-based interferometers 
have their highest sensitivity, and hdf) < 10^^^ at the lower range accessible to LISA, 
/ ~ 10"^ Hz. 



7.2 The COBE bound 



Another important constraint comes from the COBE measurement of the fluctuation of 
the temperature of the cosmic microwave background radiation (CMBR). The basic idea is 
that a strong background of GWs at very long wavelengths produces a stochastic redshift 
on the frequencies of the photons of the 2.7K radiation, and therefore a fluctuation in their 
temperature. This is known as the Sachs- Wolfe effect [224, 168 1. We first give a qualitative 
understanding of this bound, and then quote more precise numbers. 

It is conventional to expand the CMBR temperature fluctuations on the sky in spherical 
harmonics, 

6T{n) 



(120) 



1=2 m=-l 

The dipole term is dominated by the Earth motion, and therefore is omitted. The multipole 
amplitudes aim depends on the observer position r, and |a/m(r)P, averaged over all available 
positions, is the angular power spectrum. The rms fluctuations of the temperature, averaged 
over the sky, are given by 



^2/ + l 

1=2 



47r 



(121) 



where Wi depends only on the detector. COBE measures a temperature difference between 
two antennas separated by an angle 6, and 



Wi = 2[l - Pi{cose)]e 



(122) 



where Pi{cos9) are the Legendre polynomials and a is the finite width of the beam, which 
effectively cuts the resolution at ^ ~ 1/a. For the quadrupole moment, one defines 



Q 



E 

m=-2 



{\a2r. 



Air 



1/2 



(123) 



Assuming a flat Harrison-Zeldovich spectrum, i.e. spectral index n = 1 (see sect. 9.2.2), the 
COBE result is 

Q = l8±1.6nK (124) 

while fitting also n the best fit gives n = 1.2 it 0.3 and Q = 15.3^iK. In any case, from the 
quadrupole alone we get 



5T 



~ (5 - 6) X 10" 



It is not difficult to see that the COBE bound on 0,g^ has qualitatively the form 



J7gw(/) < 



Ho 
f 



gw 

6T 



(125) 



(126) 
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and holds for 

3 X 10"^^ Hz < / < 10"^^ Hz . (127) 

Let us first discuss why this is the relevant frequency range. The limitation / > 3 x 10^^^ 
Hz is just the condition that these waves are inside the horizon today. The upper limit 
instead comes out as follows. It is important to distinguish between large angle and small 
angle fluctuations. The dividing line is given by the angle ^lss such that fluctuations on 
angular scales 9 < ^lss corresponds correspond to comoving lengths that were inside the 
horizon at time of last scattering, while 6 > ^lss corresponds to super-horizon sized scales 
at time of last scattering. Small scale fluctuations are sensitive to microphysical processes 
taking place at time of last scattering; large scale fluctuations are instead insensitive to 
microphysics, because they refer to regions that were not in causal contact at that time, 
and on these angular scales the Sachs- Wolfe effect, due either to scalar perturbations or to 
tensor perturbations, i.e. GWs, is the dominant contribution to the CMBR anisotropy. 

To compute 0lsS; one first defines conformal time today, r]Q, and at the last scattering 
surface, r^LSS as 

Jo a{t) 

,LSS = p4.-=^, (128) 

Jo a{t) ^/l + ZLSS 

where we assumed a flat matter-dominated Universe today. The quantities cr/o and cr/LSS 
are the comoving size of the horizon today and at last scattering, respectively; the comoving 
distance to the last scattering surface is 

^LSS = c% ( 1 ^ ] - cr]o . (129) 

In the standard picture, where recombination occurs at Zj-^c — 1300, then zlss — 1100. In 



non-standard scenarios, zlss could be as low as O(IOO) |243]. The angle subtended by the 



horizon at last scattering is | 251 ] 



^Lss ~ ^ ~ ^= . (130) 

For a non-flat Universe, 7^ 1, this generalizes to ^lss ^ (f^oZ-^LSs)^^^- For -^lss ^ 100 
this is a good approximation to the exact formula. The deviations for smaller redshift can 
be found in |24|]. For zlss = HOO, eq. ( [l30| ) gives 6Ilss ~ 2° 



The temperature fluctuations over an angle 9 are related to aim by 

/'(|a«m|'>, (131) 



6T\ 



for 



, 200° , , 

/ ~ ^ . (132) 



Therefore multipoles up to, say, / ~ 30, are due to large angle fluctuations, and the cor- 
responding anisotropics are dominated by the Sachs- Wolfe effect, due either to scalar 
perturbations or to GWs. 



45 



The waves that were outside the horizon at the last scattering, today have a frequency, 



/<^/5Ei^^o=^10-^'Hz, (133) 

where in the last equality we have taken zlss = 0(10"^). For a non-standard reionization 
scenario with zlss ~ 100 the bound becomes approximately 3 x 10^^'' Hz. 
Let us now explain how eq. ( |126| ) comes out. A GW with a wave number 

200" \ _i 



M^) ~ [-^) Vo' (134) 
subtends an angle 6 on the last scattering surface, and creates temperature fluctuation on 



this angular scale given approximately by | 251 | 



^)^~M^LSS,Me)] (135) 

where h[zisS: ^(^)] is the characteristic amplitude that this GW had at time of last scatter- 
ing. 

We therefore must connect the value of the amplitude at time of last scattering with 
the value of /ic(/) today. One then observes that, as long as the GW is outside the horizon, 
its amplitude is constant, while, after it enters the horizon, it redshifts with the FRW scale 
factor a{t) as 1/a. This can be derived from the fact that the energy density of the GW is 
Pgw ~ f'^hcif)'^, eq. (p^); / redshifts as 1/a and, after entering the horizon, ~ 1/a^, as 
for any relativistic particle, so that hc{f ) ~ 1/a. Therefore, denoting by ahor the value of a 
at horizon crossing, 

/lc(/)today = ^^^c(/)lSS = /iLSS • (136) 

«0 OLSS ^^LSS 

During matter dominance, the scale factor depends on conformal time as a{r]) ~ 77^. Then, 
for waves that enter the horizon during matter dominance, Ohor/flLSS — (^LSs/)~^- Using 
^LSS - %/V^LSS = 2/(Fo\/^LSs), one finally gets 

/^c(/)today ^ (0)^c(/)lSS ~ {^)' {^) ■ (137) 

Using eq. ([iTl), we therefore get 



Of course, the value of {6T/T) induced by GWs cannot exceed the observed value, eq. ( 125 ), 
and this gives the bound on /igilgw 

Using the more quantitative analysis of refs. |11, ^, where the effect of multipoles with 
2 < Z < 30 is included, this bound reads 



- tr \ 2 

/io5^gw(/) < 7 X 10-11 ( ^ ) , (3 X 10-1^ Hz </< 10-16 Hz). (139) 
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This bound is stronger at the upper edge of its range of validity, / ~ SOHq ~ 10"^^ Hz, 
where it gives 

hln^if) < 10-14 , (/ ~ 10-16 Hz) . (140) 

Another issue is whether GWs can saturate this bound, or whether instead the dominant 
contribution to 5T/T comes from scalar rather than tensor perturbations; the answer is in 
general model-dependent, and we will discuss examples in sect. |9.2.2| . 

The important point that we can draw is the following: from the explicit computations 
discussed below, we will find cosmological spectra that extends from very low frequencies, 
/ ~ Hq ~ 3 X IQ-i^ Hz, up to the frequencies accessible at GW experiment and higher. 
This will be typically the case of the spectra predicted by the amplification of vacuum 
fluctuations. The frequency dependence can often be computed reliably. Then, in order for 
these spectra to be observable in the IHz-lkHz region, they must grow fast enough, so that 
they can satisfy eq. (140), and still have a sizeble value of /igOgw at higher frequencies. The 
requirement becomes even more stringent if we combine it with the nucleosynthesis bound: 
the spectra cannot keep growing until they reach the cutoff, say at 1 GHz, because otherwise 
the NS bound on the integral of /ig^^gw would be violated. Therefore it is necessary that 
there is at some stage a change of regime, so that the spectrum flattens or starts decreasing. 

Thus, it is clear that the combination of the NS bound and the COBE bound, together 
with the request of having a sizeble value of /igilgw at, say, 1 kHz, restricts considerably 
the cosmological models that can produce an interesting spectrum of vacuum fluctuations. 
Quite remarkably, there is however one candidate, theoretically well motivated, that passes 
these severe restrictions; this is the string cosmology model of Gasperini and Veneziano, 



and we will discuss it in sect. 9.3 



7.3 Pulsars as GW detectors: msec pulsars, binary pulsars and pulsar 
arrays 

Pulsars are rapidly rotating, highly magnetised neutron stars formed during the supernova 



explosion of stars with 5 to 10 solar masses (for a recent review see ref. |184| ). Soon after 
their discovery, it became clear that pulsars are excellent clocks. Actually, single pulses have 
a rather herratic profile in time, but after integrating over a number of pulses one finds an 
integrated profile that is extremely stable, and can be used as a clock. Most of the known 
pulsars have a period of the order of one second, while a growing fraction of this sample has 
periods in the range 1.5 to 30 msec. These millisecond pulsar are an extremely impressive 
source of high precision measurements. For instance, the observations of the first msec 
pulsar discovered, B1937+21, after 9 yr of data, give a period of 1.557806 468 819 794 5 ± 
0.000 000 000 000 0004 ms. The speed up of the orbital period of the double neutron star 
system B1913+16 provided experimental evidence for the existence of GWs and a test of 



General Relativity at the level of 1% |151| , 241, 242, 91]. This phenomenal stability makes 



pulsar competitive with atomic clocks (for instance PSR B1855+09 becomes competitive 
with atomic clocks after 3 yr of observations), and has even allowed to detect orbiting bodies 
smaller than the Earth from the radial acceleration that they induce on the pulsar [184]. 



Pulsars are also a natural detector of GWs ||, |2|, |T4|, |l], ||, |40|, |T6|. In fact 



a GW passing between us and the pulsar causes a fluctuation in the time of arrival of the 
pulse, proportional to the GW amplitute hc{f). Consider a GW traveling along the z axis. 
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and let 6 be the polar angle of the Earth-pulsar direction, measured from the z axis. Then 
the Doppler shift at time t induced by the GW is 



Ai/(t) 



cos I 



s s 

cos 2'ip hj^{t ) + sin 2^ (t 1 

c c 



(141) 



where, as in sect. 3, is a rotation in the plane orthogonal to the direction of propagation 
of the wave (in this case the {x,y) plane), corresponding to a choice of the axes to which 
the + and x polarization are referred. The factor (1/2) (1 — cos 9) is the standard angular 
dependence of the Doppler shift, and s is the distance along the path, so that t — s/c is the 
time when the GW crossed the Earth-pulsar direction. 

If the uncertainty in the time of arrival of the pulse is e and the total observation time is 
T, this 'detector' would be sensitive to /ic(/) ~ ^/T, for frequencies / ~ 1/T. The highest 
sensitivities can then be reached for a continuous source, as a stochastic background, after 
one or more years of integration, and therefore for / ~ 10~^ — 10~^ Hz. Based on the data 
from PSR B1855-I-09 and correcting an error from a previous analysis, ref. [^46|| gives a 
limit, at / = = 4.4 x 10"^ Hz, 



^o^gw(/*) < 1.0 X 10- 



(95%c.l.) 



(142) 



or 



^o^^gw(/*) < 4.8 X 10"^ 



(90%c.l. 



(143) 



Since the error on hc{f) is proportional to 1/T and /igilgw ~ h'^, the bound for / > is 



(at 90% c.l.) 



hln^M < 4.8 X 10-9 ( j- 



(144) 
lOV*. For 



and therefore it is quite significant (better than the NS bound) also for / 
f < f*, instead, the pulsar provides no limit at all. 

These limits are likely to continue to improve since a number of recently discovered msec 
pulsars appear to have a timing stability comparable or better than PSR B1855+09. These 



limits are especially significant for cosmic strings, as we will see in sect. (10.2). 

Another important limit comes from pulsars in binary systems, e.g. with a white dwarf 
or neutron star companion. These systems have a second clock, which is the binary orbit 
itself, and are very clean. The change in the orbital period P can be computed from 
General Relativity and in many cases is neglegibly small. This has the advantage that it 
not necessary to make a fit to an unknown derivative of the orbital period, P, and subtract 
it to isolate the contribution of GWs. This gives a limit that holds for a large range of 
frequencies, 

where T is the observation time and D the distance to the pulsar. A detailed analysis of 



this limit has been performed in |169|, correcting previous estimates. In particular the limit 



from B1855+09 turns out to be the most stringent and gives [169] 



/igl^gw < 2.7 X 10^ 



(146) 
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in the frequency range 10^^^ < / < 4.4 x 10^^ Hz. 

Finally, one can observe that these limits have been set using a pulsar as a single detector. 
From the analysis of sect. 3 it is clear that if one can use two or more pulsars as coincident 
detectors one could obtain a much better sensitivity [ [L43| , 116, 184|. From eq. ( |141| ) we see 



that the Doppler shift from the i-th pulsar has the general form 



a^sit) + n^{t), (147) 



where s{t) is the GW signal, a geometrical factor depending of the line-of-sight direction, 
and ni{t) contains the noise intrinsic to the pulsar timing. The cross-correlation between 
two pulsars then gives 

aiaj{s^{t)) + ai{sni) + aj{snj) + {niUj) , (148) 

and, since the noises from different pulsars are uncorrelated, increasing the observation 
time this tends to aiaj{s^{t)) . Furthermore, the correlation between many pulsars allows 
to subtract the errors due to terrestrial clocks, since there is now the possibility of timing 
multiple pulsars against each other. Furthermore, precise timing also require to subtract 
the effect of the Earth motion, and the precision of pulsar timing measurement is now near 
the limit of models of the Earth motion. Again, multiple pulsar timing allows to subtract 
this effect. 

Performing pulsar coincidences using the present database of long-term timing obser- 



vation does not improve on the limits obtained with single pulsars |184]. The basic reason 
is that, in order to effectively subtract the error due to the Earth motion, one needs three 
pulsars widely separated across the sky. A fourth pulsar is needed as a clock, and a fifth 



to measure the GW amplitude [116]. However, the galactic msec pulsars known around 
1990, for which we now have many years of data taking, were all concentrated in the same 
region of the sky (as a consequence of the fact that they were all found at Arecibo); in 
this situation one cannot subtract effectively the error due to the Earth motion. However, 
the msec pulsars known today are distributed much more uniformly across the sky, and 
therefore the continued timing of these pulsar in the next few years will greatly improve the 
sensitivity of pulsar arrays as GW detectors. 

Finally, we observe that the limit on h'^Q.^^ from pulsar timing holds for all sort of GW 
backgrounds, not necessarily of primordial origin, contrarily to the nucleosynthesis bound 
discussed in sect. that of course holds only for GWs that were already present at time 
of nucleosynthesis. 

7.4 Light deflection by GWs 

We finally mention that there has been quite some activity in trying to obtain bounds on 
/igJlgw from the effect that a stochastic background of GWs would have on the propagation 
of light from distant sources |265, 51]. Examining the proper motion of quasars, ref. [141| 



sets a limit 



d//i^Ogw < 0.11, (149) 
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where ~ 2 x 10^^ Hz is the inverse of the observation time, and states that considerable 
improvement in this Hmit should be possible in the next decade. Note that the limit is 
weaker than the one given by nucleosynthesis, but it is not restricted to GWs present 
before nucleosynthesis. 

A number of other mechanisms related to light deflection by GWs have been inves- 
tigated. In particular, the apparent clustering in the galaxy-galaxy correlation function 
which, according to ref. [181|, provides a bound h^Qg^ < 10^^ for wavelength between a 
few tens of kpc to a few hunderds Mpc, i.e. for 10~^^ Hz < / < 10~^^ Hz. This bound 
however has been reanalyzed in ref. | 158 | where (confirming older computations [263|) it 
has been found that, to linear order in the metric perturbation, the effect due to GWs 
does not grows linearly with the distance to the source, contrarily to previous claims, but 
only logarithmically, and therefore does not cause an appreciable galaxy clustering. Simi- 
larly, no useful bound on /igOg^ can be set from lateral displacement of caustic networks 
or weak-lensing of distant objects. This result is valid for stochastic GWs, and does not 
apply immediately to the deflection of light by GWs from localized sources. However, in 
this case the analysis of ref. [93| confirms the conclusion that these effects are too small to 
be of observational interest. 



8 Production of relic GWs: a general orientation 

Theoretical predictions of relic GW backgrounds are necessarily subject to large uncertain- 
ties, due either to the fact that we must use physics beyond the Standard Model, or to 
uncertainties in the details of the cosmological mechanisms. It is therefore important to 
understand what features of a given result are relatively general, and what are specific to a 
given model. 

In this section we discuss results on the characteristic frequency, the form of the spectrum 
and the typical intensity, that we will consider as a sort of 'benchmark'. The results from 
various specific computations can than be better understood from the comparison with these 
benchmarks, and in particular one can get a feeling of what results are relatively general 
and model independent, and what results are really peculiar to a specific computation. 

8.1 The characteristic frequency 

First of all, we discuss what is the typical frequency where we can expect to find today 
a signal of cosmological origin. Obviously, part of the answer depends on the dynamics 
of the production mechanism, and therefore is model-dependent, and part is kinematical, 
depending on the redshift from the production era. First of all, it is useful to separate the 
kinematics from the dynamics. 

To specify the kinematics, we need a cosmological model. We consider the stan- 
dard Friedmann-Robertson- Walker (FRW) cosmological model, consisting of a radiation- 
dominated (RD) phase followed by the present matter-dominated (MD) phase, and we call 
a{t) the FRW scale factor. The RD phase goes backward in time until some new regime 
sets in. This could be an inflationary epoch, e.g. at the grand unification scale, or the RD 
phase could go back in time until Planckian energies are reached and quantum gravity sets 
in, i.e., until t ~ tpi ~ 5 x 10"'^^ s. 
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A graviton produced with a frequency /*, at a time t = t^^ within the RD phase, has 
today {t = to) a redshifted frequency /o given by /o = /*a(t*)/a(io)- To compute the ratio 
a(i*)/a(to) one uses the fact that during the standard RD and MD phases the Universe 
expands adiabatically. The entropy per unit comoving volume is 

S = const. gs{T)a^{t)T^ , (150) 



where gsiT) is defined by [168| 



i=bosons i=fermions 

Here gi counts the internal degrees of freedom of the i-th particle (spin, color, etc.). 

Thus gsiT) is a measure of the effective number of degrees of freedom at temperature 
T, as far as the entropy is concerned. Another useful quantity is 

E S^(^y + l E 9^(^y (152) 
i=bosons j=fermions 

which counts the effective number of degrees of freedom at temperature T, as far as the 
energy is concerned, and that entered our discussion of nucleosynthesis, sect. |7[^. In the 
early Universe most particles had a common temperature, and g{T) and gsiT) are indis- 
tinguishable. In the Standard Model, at T>300 GeV, they become constant and have 
the value gs = 9 = 106.75, while today, assuming 3 neutrino species, gsiTo) ~ 3.91 and 



gsin) ~ 3.36 [|68| 



The sum in eqs. (151) and (|152| ) runs over relativistic species. This holds if a species 
is in thermal equilibrium at the temperature Tj. If instead it does not have a thermal 
spectrum (which in general is the case for gravitons) we can still use the above equation, 
where for this species Tj does not represent a temperature but is defined (for bosons) by 
Pi = 9ii'^'^ /'^0)Tf , where pi is the energy density of this species. Using 

gsiT,)a^iQT^ = gsiTo)a^ito)T-^ , (153) 

and To = 2.728 ± 0.002K |ll| one finds flT^] 



A.;.*)=,8,0xl0-»/.(^)"'Vi^). (154) 

aito) \gs[T*)J \ u J 

We must now ask what is the characteristic value of the frequency /* of a graviton produced 
at time t*, when the temperature was T^. Here of course the dynamics enters, but still some 
general discussion is possible. One of the relevant parameters in this estimate is certainly 
the Hubble parameter at time of production, Hit^) = H^. This comes from the fact that 
is the size of the horizon at time t^. The horizon size, physically, is the length scale 
beyond which causal microphysics cannot operate (see e.g. [168|, ch. 8.4), and therefore, for 
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causality reasons, we expect that the characteristic wavelength of gravitons or any other 
particles produced at time will be of order or smaller. | 
Therefore, we write 

A, = eH-^ . (155) 
The above argument suggests e < 1. During RD, = (87r/3)Gprad- Then 



^2 _ 87r^g,T_^ 

and, using = i/^=/e, eq. ( |154| ) can be written as | 17C ] 



Hi = " , (156) 
90M^ ' ^ ' 



1 (a^V''' 



fo ~ 1.65 X 10-^- ( ) ( — ) Hz . (157) 

This simple equation allows to understand a number of important points concerning the 
energy scales that can be probed in GW experiments. Basically, the effects of the dynamics 
have been isolated into the parameter e. 

The relation between time and temperature during the RD phase then tells us how 
far back in time are we exploring the Universe, when we observe a graviton produced at 
temperature T*, 

2.42 /MeV\2 

t.^-^i^-^j sec, (158) 

and therefore, detecting a GW that today has a frequency /o, we are looking back at the 
Universe at time 

t.:.6.6xl0--^(-^) (— ) sec. (159) 

The simplest estimate are obtained taking e = 1 in eq. (|157] ). Table 3 gives some representa- 
tive values for the production time and the temperature of the Universe T^, corresponding 
to frequencies relevant to LISA and to ground based interferometers, setting e = 1. 

However, one should keep in mind that the estimate A^, ~ -ff^T^j i-e- e ~ 1, can sometimes 
be incorrect, even as an order of magnitude estimate. Below we will illustrate this point 
with some specific examples. The values in table 3 are more correctly interpreted as a 

*On a more technical side, the deeper reason has really to do with the invariance of general relativity under 
coordinate transformations, combined with the expansion over a fixed, non-uniform, background. Consider 
for instance a scalar field (i>(x) and expand it around a given classical configuration, (\)(x) = (/)o(a;) + &^{x). 
Under a general coordinate transformation x — > a;', by definition a scalar field transforms as (\>{x] ^'(a;') = 
<^{x). However, when we expand around a given background, we keep its functional form fixed and therefore 
under a; — > x', 00 (a;) — > 0o(a;')i which for a non-constant field configuration, is different from 00 (a^)- It follows 
that the perturbation S(ji{x) is not a scalar under general coordinate transformations, even if 4i{x) was a 
scalar. The effect becomes important for the Fourier components of 5(j){x) with a wavelength comparable 
or greater than the variation scale of the background <l)o{x). (We are discussing a scalar field for notational 
simplicity, but of course the same holds for the metric tensor gfyj^). In a homogeneous FRW background 
the only variation is temporal, and its timescale is given by the . Therefore modes with wavelength 
greater than are in general plagued by gauge artefacts. This problem manifests itself, for instance, 

when computing density fluctuations in the early Universe. In t his case one flnds s purio us modes which can 



be removed with an appropriate gauge choice, see e.g. ref. [168], sect. 9.3.6 or ref. [202| 
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detection frequency 


production time 


production temperature 


10"'' Hz 


7 X 10"^^ sec 


600 GeV 


1 Hz 


7 X 10"^^ sec 


6 X 10** GeV 


100 Hz 


7 X 10"^^ sec 


6 X 10« GeV 


10^ Hz 


7 X IQ-'^' sec 


6 X 10^ GeV 



Table 3: The production time and the production temperature T* for GWs observed 
today at frequency /o, if at time of production they had a wavelength of the order of the 
horizon length. 



starting point for understanding how the result is affected by the dynamics of the specific 
production mechanism. 

An important question is whether GW experiments can explore the very high energy 
regime, when the Universe had temperatures of the order of a typical grand unification 
scale, or of the order of the Planck scale. 

The scale of quantum gravity is given by the Planck mass, related to Newton constant by 
G = 1/Mpy More precisely, since in the gravitational action enters the combination SttG, 
we expect that the relevant scale is the reduced Planck mass 

Mpi/(87r)i/2 ^ 2.44 x 10^^ 

GeV. Using eq. ([TSTD with = Mpi/(87r)i/2 and e = 1 gives 



/o~400 f^j GHz, (T, = Mpi/v^). (160) 



The dependence on g'^, is rather weak because of the power 1/6 in eq. ( |157| ). For = 1000, 
/o increases by a factor ~ 1.5 relative to gi* = 100. Instead, for = Mqut ~ 10^^ GeV, 
and 5* ~ 220 (likely values for a super symmetric unification), 

/o~2GHz, (T, = Mgut). (161) 

Using instead the typical scale of string theory, Ms, the characteristic frequency is between 
these two values, since Ms is expected to be approximately in the range W^^ <,Ms / Mpi<,10~^ 



|161]. 



Thus, the expected peak frequency of GW spectra produced at these very early epochs 
is in the GHz region. In this region no present or future experiment can operate, and the 
reason is easily seen from eq. (^2|): considering that the maximum value of /lo^^gw is fixed 
by nucleosynthesis, the required sensitivity in the GW amplitude goes as f~^^'^. At / = 
1 GHz one needs a value oi hj smaller by a factor 10^, compared to a value that would give 
the same sensitivity in /ig^gw at / = IkHz. 

However, one should not hurry toward pessimistic conclusions. The characteristic fre- 
quency that we have discussed is the value of the cutoff frequency in the graviton spectrum. 
Above this frequency, the spectrum decreases exponentially, and no signal can be detected. 
Below this frequency, however, the form of the spectrum is not fixed by general arguments. 
We discuss the issue in the next section. 
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8.2 The form of the spectrum 

For a thermal spectrum, the occupation number per ceh of the phase space, n/, is given by 
the Bose-Einstein distribition and therefore, below the characteristic frequency, 

1 kT , , 

= ,2.f/kT _ 1 - ^ (162) 



Then, from eq. (26) we then find that at low frequencies a thermal spectrum corresponds 
to 

17g„(/)~/3. (163) 

Of course, if we have a thermal spectrum with a cutoff frequency in the GHz region, the 
value of /igr2gw(/) at, say, / = 1 kHz is utterly neglegible. However, we have seen in the 
Introduction that below the Planck scale gravitons interact too weakly to thermalize, and 
therefore there is no a priori reason for a ~ /'^ dependence. The gravitons will retain 
the form of the spectrum that they had at time of production, and this is a very model 
dependent feature. However, from a number of explicit examples and general arguments 
that we will discuss below, we learn that spectra flat or almost flat over a large range of 
frequencies seem to be not at all unusual in cosmology. In particular: 

1. The amplification of vacuum fluctuations always give spectra that extend over a huge 
range of frequency, from 3 x 10^^^ Hz, corresponding to wavelength equal to the 
present size of the horizon, up to a cutoff in the GHz region. The spectrum found in 
some cases in string cosmology is particularly interesting from this point of view, since 
it can satisfy the COBE and nucleosysnthesis bound, and still have a large intensity 
at the frequency accessible to experiments. The basic reason for this behaviour is the 
fact that in the problem there is only one relevant lengthscale, given by the value of 
the Hubble parameter. 

2. Another example of relic GW spectrum which is almost flat over a very large range of 
frequencies is provided by GWs produced by the decay of cosmic strings, discussed in 
sect. |ia^. In this case the spectrum has a peak around / ~ 10 Hz, where /ig^gw can 
reach a few times 10~^, and then it is almost flat from / 10~* Hz to the GHz region, 
where it has the cutoff (fixed by the arguments previously discussed, considering that 
the relevant scale for cosmic string is Mqut)- The basic reason for such a behavior 
is the scaling property of the string network, which says that a single lengthscale, 
the Hubble length, characterizes all properties of the string network |254, ^. The 



network of strings evolves toward a self-similar configuration, with small loops being 
chooped off very long strings, and the typical radiation emitted by a single loop has 
a wavelength related to the length of the loop. 

3. Phase transitions can also give spectra that extend over a large range of frequencies; 
in this case, as we will see in sect. (|10.1.3|) , the reason is that in an expanding Universe 



a field cannot have a correlation length larger than the horizon size, basically because 
of causality. 

These facts have potentially important consequencies. They suggest that, even if a 
spectrum of gravitons produced during the Planck era has a cutoff at frequencies much larger 
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than the range of frequencies accessible to interferometers, still in this range we can hope 
to observe the low-frequency part of these spectra. 



8.3 Characteristic intensity 



The numbers discussed in sect. ( [7.1D give a first idea of what can be considered an interesting 
detection level for /io^gw(/)j which should be at least 



/io^^gw ^ 5 X 10"'^ , (164) 



especially considering that the bound (118) refers not only to gravitational waves, but to all 
possible sources of energy which have not been included, like particles beyond the standard 
model, primordial black holes, etc. 

The next question is whether it is reasonable to expect that some cosmological pro- 
duction mechanism saturates the nucleosynthesis bound. This of course depends on the 
production mechanism, but some relatively general considerations are possible. 

First of all, it is clear that, if GWs are produced at the Planck scale by collisions and 
decays together with the photons that we observe today in the CMBR, and there is not an 
inflationary phase at later time, we expect roughly ~ p^. Since hlp^/pc ^ 2.481 X 10"^, 
in this case the bound ( |115| ) is approximately saturated. More precisely, if at some time 
t = both the photons that we observe in the CMBR and gravitons were produced, and 
with a total fraction of the energy density in gravity waves O*, then, taking into account 
the different redshifts of GWs and photons due to the fact that GWs decouple immediately, 
at the present time we would have |170f| 



/igOgw ~ 1.67 X 10"^ (—] 



1/3 



\ g* J 



. (165) 



If instead the mechanism that produces GWs is different from the mechanism that produces 
the photons in the CMBR, often it is still possible to relate the respective energy densities, 
simply because the scales of the two processes are related. A few examples will illustrate 
the point. 

1. In sect. we will discuss the spectrum produced in string cosmology by the am- 
plification of vacuum fluctuations. However, the typical value of the characteristic 
frequency and the peak value of /igOgw can be understood using only general argu- 
ments, following ref. ||57|| . 

The production temperature is fixed by the criterium H^, ~ Ms, since in this model 
the Hubble parameter is stabilized by stringy corrections, and therefore at the string 
scale, see sect. |9.3i Using if* ~ T^/Mpi, this corresponds to an effective 'temperature' 
at time of production ~ (MsMpi)^/'^ = (M5 /Mpi)V2Mpi. Redshifting T* we get a 



characteristic frequency today 
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and a peak energy density 



dlogf \MpJ dlog/ ^ ' 

Therefore is related to p^, although with a suppression factor (Ms/Mpi)^ which 
is expected to range between 10~^ and 10~^, and numerical factors. 

2. Another mechanism that we will discuss is the production of GWs through bubble 
collisions when inflation terminates with a first order phase transition. In this case 
there is only one energy scale, the vacuum energy density M during inflation (so 
that is the false- vacuum energy density). The value of M fixes the reheating 
temperature, and therefore p-y, from {TT'^g^/30)T^Yi — -^^) ^-^d fixes also the energy 
liberated in gravitational waves. Therefore the energy density in photons and in 



gravitons are related. Writing, as in sect. 3.1, A^, = ei^^ , where A*,-ff* are the 



typical wavelength produced and the Hubble parameter at time of production, the 



computation of sect. 10.1.1 shows that, in a strongly first order phase transition, 

Pgw ~ e^P-y , (168) 
so that again and p-y are related. 

From these explicit examples we see that, independently of the production mechanism, 
an equation of the form eq. (|168| ) is quite general. The scale for the intensity of a relic 
GW background is indeed fixed in many cases by p^, which therefore gives a first order 
of magnitude estimate of the effect. However, there are also suppression factors, like the 
factor in eq. ( 237 ) or (Mg/Mpi)^ in eq. ( [1671 ), that, together with the exact numerical 
coefficients, are crucial for the detection at present experiments. A value e ~ 0.1 would 
allow detection at the level /ig^^gw ~ 10~^ while e ~ 10~^ would require h^Qg^ ~ 10~^, 
which is beyond the possibilities of first generation experiments. 

Furthermore, another suppression factor is present if we have a spectrum of relic GWs 
with a total energy density pg„ ~ p^, and a cutoff in the GHz, and we want to observe it in 
the kHz region. In this case, as we discussed, we must hope that the spectrum is practically 
flat between the kHz and the GHz. While we have seen that there are various examples 
of such a behaviour, the fact that the energy density is spread over such a wide interval of 
frequencies diminuishes the maximum allowed value of hQilgvj{f) at a given frequency, since 
the nucleosynthesis bound is a limit on the integral of hQQg^{f) over all frequencies. The 
nucleosynthesis bound (|118|) then gives a maximum value at IkHz 



hlng^iimz) < 5-6xlO-6(Ar 3) ^ ^ ^ _ ^^^^^ 

If the spectrum extends further toward lower frequencies, the maximum value of h^Qg^f 
decreases accordingly, with a factor log (1 GHz //mm) instead of log(lGHz/lkHz). 
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9 Amplification of vacuum fiuctuations 



9.1 The computation of Bogoliubov coefficients 

The amplification of vacuum fluctuations is a very general mechanisms, first discussed in a 
cosmological setting in refs. [ [L38| , |236| ] and then in many other papers, see e.g. [ |223| , 106, 
1) 0! 0) HI) 114, 139 1 . Specializing immediately to a FRW metric, it is convenient to 



introduce conformal time ry, related to cosming time thy drj = dt/a{t), so that 

ds"^ = a^{rj) (-dif + dx^) . (170) 

(In the following we set c = 1). In the presence of a classical GW in a FRW background, 
we write 

g^Lu = a^iv) iVfiu + hf,u) , (171) 
with r]fj_iy = (-, +, +, +) and, in the TT gauge. 



A=+,x k 

Note that x are the comoving coordinates; the physical coordinates are Xp^ys = a(t)x. 
Correspondingly, k is a comoving momentum, and the physical momentum is kp^yg — 
k/a{t). 

The Einstein equations, linearized over the FRW background, give an equation for <p^, 

< + 2^(A;, + A;2^k = 0, (173) 

where the prime is the derivative with respect to r], and we omitted the index A. This 
equation is just the Klein-Gordon equation for a homogeneous scalar field in a FRW back- 
ground, 

□0 = -^d^ iV^g^^d,) = 0. (174) 



The factor ^/SttGn in eq. ( 172 ) is a convenient choice of normalization, such that the action 
of (f)^ derived from the linearization of the Einstein-Hilbert action corresponds to the action 
of two (^ = -|-, x) canonically normalized scalar fields. Therefore one can use the standard 
technique of Bogoliubov coefficients for a scalar field |44]. One considers the situation in 
which, over a timescale AT, there is a change of regime in the cosmological evolution, 
between two regimes that we denote I and II (for instance from an inflationary regime to a 
standard radiation dominated era). Let be the time of transition between the two regimes, 
and consider a mode whose physical frequency, at time t*, is /*. For values of /* such that 
27r/*AT <C 1, the change of regime can be considered abrupt, while if 27r/*AT ^ 1 the 
change of regime is adiabatic. It is convenient to define a new variable 

Mv) = -Mv) ■ (175) 
a 

In terms of (f>\i, eq. ( |173|) becomes 

V'k+ ^'--Vk = 0. (176) 
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Let us denote by /k(??) the solution of this equation with the scale factor of phase I, and 
-^k(^) for phase II. Then the mode expansion of hat in the phase I is 



Kb = V^ttGn 



1 



(27r)3V2fc a{ri) 



aA{k)Uv)e^^-- + a\ik)f^i^)e-^^-- e^.^ (177) 



(we have used the fact that our e^j, are real, eq. (P)); a^(k),a|4(k) are the creation and 
annihilation operators in this phase, so that the vacuum in phase I is defined by 



a+(k)|0)7 = ax(k)|0)7 = 0, 



(178) 



and we can construct the Fock space appropriate to this phase acting with a) (k) . Similarly, 
in phase II, one has the expansion 



ab 



1 



and a new vacuum state \0)ii such that 

A+{k)\0)ji = A^{k)\0)ji = 0. 



^^(k)Fk(r?)e^''- + 4(k)Fk*(r?)e-H e^,(A) (179) 



(180) 



Since both the (/k, f^) and the (Fk, F^) are a complete set, we can express one in terms of the 
others; of course, since the /k alone are not a complete set, the expansion of -Fk will involve 
both /k and and therefore there will be a mixing of positive and negative frequency 
modes. The relation between Fk and (/k, /k) is known as a Bogoliubov transformation, 



Fk = Yl («kk'/k' + /3kk'/k') • 
k' 



(181) 



Inserting this relation into eq. (179) and comparing with eq. (177) one finds the relation 
between creation and annihilation operators (we omit hereafter the index A = +, x), 



and 



a(k) = ^[ak'kA(k') + /3k'k^^(k') 
k' 



The Bogoliubov coefficients satisfy 

["kikOkak - /5kik/3k2k] 



(^kiko 



Y ["kik/3k2k - /^kikOkak] = 



(182) 

(183) 

(184) 

(185) 



If the cosmological background metric is time-dependent but isotropic and spatially homo- 
geneous, the gravitational field can give energy to the particles, but not momentum, so 
that 



Okk' 
/3kk' 



a/(5kk' 
/?/5kk' • 



(186) 
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Eqs. ([18^JT83|) give the formal relation between the Fock spaces constructed in phase I and 



II. Now, the basic physical idea is the following: suppose that just before the transition the 
quantum state was |s). For instance, we can express \s) in terms of the occupation numbers 
Uf, relative to the creation operators cl^jClj- Then, for modes such that the transition is 
sudden, that is 27r/*AT <C 1, the physical state does not have time to change during the 
transition, so that it is still given by \s) = just after the transition. However, we now 

must express the occupation numbers with respect to Nf = J^^Af. Using the Bogoliubov 
transformation, one immediately finds 



Nf=nf + 2\(3f\'[nf + -^ . (187) 

This shows, first of all, that any preexisting value of n/ is amplified, with an amplification 
factor 1 + 2|/3jp. And furthermore, even the 'half-quantum' due to vacuum fluctuations is 
amplified: due to the mixing between positive and negative frequencies, even the vacuum 
state of phase I is a multiparticle state in the Fock space appropriate to phase II. 

The above discussion refers to frequencies /* such that 27r/*AT ^1. In the opposite 
limit 27r/*AT » 1, the mode whose physical frequency, at t = i*, is /*, sees the transition 
between the two regimes as adiabatic, and the quantum state has the time to follow smoothly 
the evolution of the scale factor. In this case, therefore, there is no particle production, or 



more precisely the amplification factor approaches one exponentially [44|. Therefore, for 
any practical purpose, there is a cutoff in the spectrum produced by the amplification of 
vacuum fluctuations, at 

/r^" ^ . (188) 

To obtain the value of the cutoff frequency today, one then redshifts /™^^ from time of 
production t = t^,to the present time to- Iii a cosmological setting, the scale of time variation 
is given by the Hubble constant, so that AT ~ H~^; then the condition 27r/^,Ar » 1 is 
met when the reduced physical wavelength A*/(27r) is smaller than the horizon; thus, modes 
that at time of transition where outside the horizon are amplified. Instead, no amplification 
occurs for sub-horizon sized modes. 

It is also interesting to note, from eq. (187), that even if we have an inflationary phase. 



the occupation numbers Nf after inflation are very sensitive to the occupation numbers nj 
before inflation. This is in contrast with what happens in most cases after an inflationary 
phase: usually inflation practically erases the information about a previous era. This does 
not happen for the amplification of vacuum fluctuations, and the technical reason behind 
this is that the quantity which is amplified is a number of particle per unit cell of the phase 
space. The volume of a cell of the phase space, xd?k / {2t:)^ , is unaffected by the expansion 
of the Universe, contrarily to a physical spatial volume d?x ~ a?{t). 



9.2 Amplification of vacuum fluctuations in inflationary models 

We now compute the Bogoliubov coefficients for the transition between an inflationary 
phase and the radiation-dominated (RD) phase, followed by the matter-dominated (MD) 
phase. We start from the simplest case of exact De Sitter inflation and then we compute 
the modifications for slow-roll inflation. 
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The amplification of vacuum flucutations at the inflation-RD transition is the most 
studied example of this general phenomenon |l38|, |236|, |22|, |ll^, |, [I], |ll|, |2|, |3^, |l40|| , 



but, as we will see, it is by no means the most promising from the point of view of the 
detection at interferometers or resonant masses. Still, we start with this example because 
is the simpler setting to understand the computational technique. 

9.2.1 De Sitter inflation 

We start with a FRW cosmological model with ds^ = dt^ — a?{t)d'x^ ^ and we concentrate for 
the moment on the particle production at the De Sitter-RD transition. The scale factors 
are a{t) ~ e^* for — oo < t < ti, and a{t) ~ t^^"^ for ti <t < tcq, where and teq is the time 
when MD sets in. In terms of conformal time t], 

ain) = (189) 



for — cxD < r] < rji, with r]i < 0, and 



1 



<V) = -ET^iv - 2m) , (190) 



for rji < rj < r]cq- We have fixed the constants so that a{r]) and a'(r/) are continuous 
across the transition. Here we are considering an istantaneous transition at i = ti. More 
realistically, we should consider a transition that takes place in a time At. As already 
discussed in the previous section, as long as we consider frequencies such that 27r/*At <^ 1, 
the approximation of an istantaneous transition gives neglegible errors. When 27r/*At ~ 1, 
the approximation breaks down, and there is an exponential suppression of the particle 
production. 

The solution of eq. (|176| ) with these scale factors is easily found. In general the solution 
is a superposition of a term that oscillates as e~^^^ and a term e^^"^. For < r/i, we impose 
the boundary condition that only the positive energy solution, i.e. the term ~ e~^^^, is 
present. This is not an assumption on the state of the system at ij < rji, but simply a trick 
to compute the Bogoliubov coefficient, which can be defined by matching a solution with 
only positive frequencies at r] < rji with the most general solution at r] > rji. Then, the 
solutions are 

for — oo < r] < rji, with rji < 0, and 

Mv) = {ake-'"^ + Pke''"^) (192) 
for r]i < 7] < rjcq. Requiring that the 0^, (p'k are continuous across the transition gives 

''' = ^-fk~^i' = ^^^^^ 

Note that |afcp— |/?A;P = 1- The number of particles produced at the De Sitter-RD transition, 
per cell of the phase space, is given by eq. (187); in particular, if before the transition the 
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system was in the vacuum state, = 0, after the transition 

As discussed before, k is the comoving momentum and is related to the physical momentum 
fcphys by /cphys = k/a{t). We must now express in terms of physical quantities observed 
today, taking into account the redshift during the RD and MD phases. Of course, during 
the RD-MD transition there will be a further creation of quanta; we will discuss it below, 
and for the moment we express the spectrum of particles produced at the De Sitter-RD 
transition in terms of physical quantities measured today. 

We reserve the notation / for the physical frequency observed today, so that 

27Tf = (A:phys)o = , (195) 

where to is the present value of cosmic time. Therefore, 

fcM=2.Mto)M = ^|^ = ^(f^) (^) , (196) 




where in the second equality we have used a{ti) = 

(—\ EE 1 + zeq = 2.32 X 10^ (l^o/io)^2T75 - 2-40 X 10^ (i^ohl) , (197) 

\ icq / 

is the total density of the Universe in units of the critical density, and we have taken 
T = 2.728 K. For the time of matter-radiation equilibrium we have [168| teq = 4.3608 x 
10iO(17o/ig)-2(T2.75)^ ^ 4.1 • IQio n^^h^^s. 

From eq. (|196| ) we see that the parameter rji can be traded for a more meaningful 
parameter /i defined by 

k\vi\ = ^, (198) 



so that 



Note that here ti is not a free parameter: during RD, H{t) = l/(2t), so that H{ti) = 
l/(2ti), since ti is the time when RD begin; and since the scale factor and its derivative are 
continuous across the transition, H{ti) is equal to the constant value H during the De Sitter 
phase, so that ti = \/{2H). Then, inserting the numerical values, 

/i^lO'f -,n-f.. V^' Hz. (200) 



,10-4Mpi, 

Note that both /iq and JIq cancel in /i. We have chosen lO^^Mpi as a useful reference value 
for for reasons that will be clear below. In terms of the physical frequency observed 
today, /, and of the parameter /i, 

^f=\{jy^ (201) 
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if before the transition n/ = 0. Inserting this value into eq. ( Pq ) we see that / cancels and 
we get a flat spectrum, 

^o^^-(^)-10""(t^)'- (202) 

This is the contribution of the De Sitter-RD transition. A further amplification comes from 
the RD-MD transition. The mode that, at the time of the transition between RD and MD, 
had a physical frequency much bigger than the typical timescale of the transition, are not 
subject to further amplification. The typical timescale of the transition is given by the 
inverse of the Hubble constant at time of equilibrium, which is H{teq} = 4:{^/2 — l)/(3teq) 
(see [168], pg. 60). A mode that, at t = teq, had a physical momentum A^phys = H{teq), 
today has a physical frequency /eq given by 

/eq = ^^4^ — " ^^-''M) Hz . (203) 

ZTT O^cq ^cq 

The modes that today have a frequency / > /eq have not been further amplified at the 
RD-MD transition. Another useful way to express this result is the statement that modes 
that, at the time of transition, had a physical wavelength much smaller than the horizon 
are not amplified. Therefore, the spectrum at / > /eq is given by eq. ( |202| ). 

For 3 X 10^^^ Hz < / < /eq Hz, instead, we must compute the Bogoliubov coefficients 
for the RD-MD transition, and the total Bogoliubov transformation is obtained with the 
composition of the De Sitter-RD and RD-MD transformations. As always, the condition 
/ > 3 X 10~^^ Hz comes from requiring that the mode is inside the horizon today. Thus, 
these are the modes that were still outside the horizon at the RD-MD transition and re- 
entered the horizon before the present time. 

A straightforward computation, similar to the one performed above, shows that in this 
frequency range the spectrum now goes like 1//^ rather than being flat. So the final result 

if 3 X 10-18 Hz < / < /eq ~ lO-i*^ Hz, and 

hl^.M^lO-^'{j^^\ (205) 



for /eq < / < /i- Finally, from the definition of /i, eq. (198), we see that a mode that today 
has / > /i, at the time of transition between De Sitter and RD had a wavelength smaller 
than the horizon; therefore these modes are never amplified, so that /i is the upper cutoff 



of the spectrum. The spectrum is shown in fig. |12 . 

Note that the behaviour 1//^ in the region 3 x 10^^^ Hz < / < 10^^^ Hz is the same 
as the COBE bound, and is in fact due to the same physical reason, i.e. to the further 
amplification of modes that are inside the horizon at time of matter-radiation equilibrium. 
The comparison of eq. (204) with the COBE bound, eq. (139), then gives an upper limit on 



the Hubble constant during De Sitter inflation |235| ] ; from a precise statistical analysis, the 
limit at 95% c.l. is 

H <6x lO^^Mpi . (206) 
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Figure 12: The spectrum of amplification of vacuum fluctuations produced by a phase of 
De Sitter inflation (solid line), with a value of H that saturates the COBE bound. The 
nucleosynthesis bound (dotted line) and the pulsar bound (triangle shaped) of fig. 11 are 
also shown for comparison. 
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9.2.2 Slow-roll inflation 



In the above section we have considered a simplified model with an exact De Sitter inflation. 
If inflation is driven by a scalar fleld slowly rolling in a potential V{(p), the Hubble parameter 
during inflation is not exactly constant, and this results in a small tilt in the spectrum. 
Furthermore, scalar perturbations will also contribute to the COBE anisotropy, and the 
ratio of tensor to scalar perturbation is fixed by the choice of V{4>)- Then, in the region 
/eq < / < A) one finds that /io^^gw(/) is not flat, as in eq. ( ^051) , but has a frequency 
dependence 

,(/) ~ n , (207) 



where is called the spectral index for tensor perturbations. If 14 is the value of the 
inflationary potential when the present horizon scale {k = Hq) crossed the horizon during 
inflation, and is the first derivative of the potential at that point, then, to lowest order in 



the deviation from scale invariance, and for single field models with a smooth potential [178, 
249| , [18^ 

„, = _^(YL)\ ,208) 



Svr VK 



Therefore nj- < 0, and the spectrum is decreasing rather than flat. Since the slow-roll 
condition is just ( v^j I ^ then \nT\ ^ 1- It is convenient to define the scalar (S) 
and tensor (T) contributions to the quadrupole anisotropy 



S 
T 



47r 
5(|aLP) 



2.2(K/M^i) 
(MpiKVV;)2 

0.61(K/M|i) 



(209) 
(210) 



where the second equalities hold for slow-roll inflation. From eqs. ( 208[ ), ( p09| ) and (210) it 
follows that 

»r = 4|. (211) 

To obtain the precise frequency dependence, we must also take into account the fact that 
the behaviours /ig^^gw ~ 1//^ and Hq^, gw ~ const, in eqs. ( |204D and ( p05D hold only in the 
limits / <C /eq and / ^ /eq. Performing accurately the matching at / = /eq, taking into 
account the tilt nx and defining keq = '2^7^ feci, finds the spectrum [124^ 



PI) 



dlogk {k/Hc, 



\2~nT 



4 k 5 k 
^+3A- + 2(a- 



eq 



(212) 



Defining = {87:/3){V^/Mp^) and inserting the numerical values, we can rewrite it as 



0.8 X 10 



-13 



lO-^Mpi 



2 /5 X IQ-^'^hnRz 



— Tlx 



^ ^ /5 X 10-i7(/i§J]o) Hz\ /6 X 10-i7(/igOo) Hz' 



(213) 
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The relation between ut, T and S, eq. (211), goes in the wrong direction for the detection at 
interferometers: the higher the contribution of tensor perturbations at COBE, and therefore 
the maximum allowed value of the GW spectrum at, say, / = 10~^^ Hz, the steeper is the 
decrease in the spectrum with the frequency. Furthermore, when extrapolating through a 
very large range of frequencies, one must also take into account that nx itself is not constant, 



but rather [172 



"-r .„^Mk (Yl) . (214) 



dlogA; Air 

Combining these effects, one finds that, at a fixed frequency / = 10~^ Hz, relevant for 
LISA, the spectrum is maximized for n-p ~ —0.025 |249| ] , or T/S" ~ 0.175. Values of this 
order are realized in various models of inflation, although there are as well models that 
predict T/S ~ 10~'^. As discussed in sect, [t!^, GWs can contribute to the anisotropy only 
on angular scales larger than approximately 1°. Combining the results of COBE with other 
observations on smaller angular scales, one can place upper limits on T/S, and the analysis 
of ref. [|9|, m gives T/S < 0.5, at 95% c.l. 

Even with the value of ut that maximizes the signal at / = 10~^ Hz, the expected 
signal at LISA is of order 10~^^, therefore various order of magnitudes below the expected 
sensitivity. 



9.3 Pre-big-bang cosmology 
9.3.1 The model 

In recent years, a cosmological model derived from the low-energy effective action of string 
theory has been proposed |253| , |13d|| (see ref. [179] for a recent review). The simplest starting 
point is the effective action in the metric-dilaton sector; at lowest order in the derivatives 
and in e"^ it is given by 



5. 



eflf 



-9 



{R + d^ct)d''(t>) 



(215) 



where (p is the dilaton field. This action receives two kinds of perturbative corrections: 
corrections parametrized by e*^, which are just higher loops in the field theory sense. And 
corrections proportional to a', where = \/2a' is the string length, which are genuinely 
string effects, and take into account the finite size of the string. Since a' has dimension 
of (length)^, it is associated with higher derivative operators as for instance the square 
of the Riemann tensor, Rf^up^- The remarkable fact observed in refs. p53| , |130[ is that 
there exits a regime in the early Universe evolution of this model, where the perturbative 
approach is well justified, and in a first approximation one can use eq. ( pl5D , neglecting the 
corrections. This regime occurs if we take as initial condition a Universe at weak coupling 
and low curvatures. Then, a generic inhomogeneous string vacuum shows a gravitational 
instability [S2[. Specializing for simplicity to a homogeneous model, the solution of the 



equations of motion derived from the action (215) reads, for generic anisotropic FRW scale 
factors ttiit), i = 1, . . . D — 1, 



Mt) = 00 + Co log(-t) , 



(216) 
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with 

D-l D-l 

5] c? = 1 , J2ci = l + co. (217) 

i=l 1=1 

This is just a generahzation in the presence of the dilaton of the Kasner solution of general 
relativity. Considering for simplicity an isotropic model, one has in particular a solution with 
Ci = —l/\/D — I, Co = —1 — \/D — I. This corresponds to a superinflationary evolution: 
the Hubble parameter H = a/a grows until it reaches the string scale ~ l/A^, and finally 
formally diverges at a singularity. When H ~ 1/As, of course, the description based on 
the lowest order action is not anymore adequate. Much work has been devoted recently 
to understanding how string theory cures the apparent singularity, and if it is possible to 
match this 'pre-big-bang' phase to standard post-big-bang cosmology |60, |160| , |101| , |88| , 



128, 187, 114, 115]. Apart from the phenomenological aspects concerning GWs that we 



will discuss, it is clear that the problem has great conceptual interest, since it a theoretically 
well motivated attempt to cure the Big-Bang singularity. 

One can then investigate whether pertubative corrections succeed in turning the regime 
of pre-big-bang accelerated expansion into the standard decelerated expansion. Indeed, the 
inclusion of a' corrections can turn the unbounded growth of the curvature into a De Sitter 
phase with linearly growing dilaton 128 |. Since in this new regime the dilaton keeps growing 
linearly, finally also e'^ gets large, even if the evolution started at very weak coupling; at 
this stage loop corrections are important, and they can trigger the exit from the De Sitter 



phase py, 115|. An example of this behaviour, using loop corrections derived from orbifold 



compactifications of string theory, has been found in ref. [115|, and is shown in fig. |l3|. 

These results are encouraging, but various technical limitations in the analysis have also 
indicated that most probably the true resolution of the singularity problem lies outside 
the domain of validity of the perturbative approach |11£]. In this case one must take into 
account that at strong coupling and large curvature new light states appear and then the 
approach based on the effective supergravity action plus string corrections breaks down. The 
light modes are now different and one must turn to a new effective action, written in terms 
of the new relevant degrees of freedom. In particular, at strong coupling D-branes [ plO|| 
are expected to play an important role, since their mass scales like the inverse of the string 
coupling, ~ 1/(7 and they are copiously produced by gravitational fields [19C|. 

For the purpose of computing the amplification of vacuum fluctuations in this model, 
the approach that can be taken is to consider a model with three phases: first an isotropic 
superinflationary evolution described by eq. (216), with the value of the constants Ci,co 
determined by eq. (|217] ). Then this phase is matched to a De Sitter phase with linearly 
growing dilaton, which may be considered as representative of a typical solution in the large 
curvature regime; the constant values of H, (j) in this phase are taken as free parameters, of 
the order of the string scale. Finally, we match this phase to the standard RD era. Of course, 
this model can only be taken as indicative of possible behaviours in string cosmology, since 
what happens in the large curvature phase is certainly much more complicated. However, 
we will see that some important characteristics of the GW spectrum depend only on the 
low curvature regime, where we can use the lowest order effective action and therefore is 
under good theoretical control. 

Variants of this model can of course be constructed. In particular, ref. | 197 | considers 
the effect of a second burst of inflation, or of the possibility that the Universe becomes 
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temporarily dominated by long-lived massive particles, as for instance the moduli fields of 
string theory, which then decay, restoring a standard RD phase 137 |. 



9.3.2 Production of GWs 



We can now discuss the production of GWs through amplification of vacuum fiuctuations 

m, |io|, |i88|, [i9|, m, m, m, n, iisi ). 



in this model |131, |126| , 56 1 (see also refs. [|^, 57 



We specialize for illustration to a isotropic model compactified to four dimensions. Writing 
the scale factors and the dilaton in terms of conformal time, from eq. ( |216| ), we find that 
for —oo<r]<ris (with r/^ < 0), we have a dilaton-dominated regime with 



1 



r/ - (1 - a)ris 



s - 7 lo; 



ar]s 
(1- 



o:)ris 



(218) 
(219) 



and the Kasner conditions eq. ( |217D give a = 1/(1 + \/3), 7 = \/3. 

At a value i] = rjg the curvature becomes of the order of the string scale, and we assume 
a De Sitter expansion with linearly growing dilaton; in terms of conformal time, this means 

1 



a(r?) 



2/? log 



7] 



(220) 



The parameter Hg is of the order of the string scale, while /? is a free parameter of the 
model. Note that 2/3 = (p/Hg, where is the derivative with respect to cosmic time. 

The stringy phase lasts for rjs < rj < rji (again rji < 0), and at r/i we match this phase 
to a standard RD era. This gives, for rji < t] < r]r, (with r]r > 0) 



a{r]) 



iv - 27?i 



(221) 



After that, the standard matter dominated era takes place. We have chosen the additive 
and multiplicative constants in a{rj) in such a way that a{rj) and da/dtj (and therefore also 
da/dt) are continuous across the transitions. 

The equation for the Fourier modes of metric tensor perturbations for the two physical 
polarizations in the transverse traceless gauge is [126| 

'^''^' + U2-y(r?)l Vfc = 0, (222) 



drf 



V{ri) 



1 



d^ 



... , . (223) 
a drj'^ V / 

Note that, for constant V{r]) reduces to eq. (|176D . Inserting the expressions (I2I8N221I ), 
the potential is 

1 



V[r,) = \ (4/.2 
1/(77) = 0, 



1) - (1 - 



-00 < T] < rjs 

r]s Kr] <r]i 
m <V <Vr 



(224) 
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where 2// = |2/3 — 3|,2i^ = |2a — 7 + 1|. From a = 1/(1 + \/3) and 7 = \/3, we get u = 0. 
The exact solutions of eq. ( ^221) in the three regions are 



r?-(l-a)r?,|C7i7f (A:|r?-(l-a)r/,|), 



Si e 



ifcr; 



B_ e 



-00 < 7] < rjs 

r]s <v <'ni 



(225) 



where hI^'"^^ are Hankel's functions. As we have seen in sect. p. 2.1 , it is convenient to trade 
i]s,rji for more physically meaningful quantities. As before, we then write 



k\r]i\ = 27r/a(to)|??i| 



27r/ a(to) ^ 27r/ 

Hs a{tl) Hs \tcq 



2/3 



t \ 
''eq 



(226) 



and therefore the parameter r/i can be traded for a parameter /i defined by 



k\rji\ 



f_ 

h 



fi ~ 4.3 • lO^^Hz 



0.15 Mpi 



1/2 



(227) 



where we have chosen a reference value Hg of the order of the string scale, and ti ~ A^. 
Similarly we can introduce a parameter fg instead of r/s, from k\ris\ = f / fs- This parameter 
depends on the duration of the string phase, and therefore, contrarily to /i, is totally 
unknown, even as an order of magnitude. However, since < \ris\, we have fs < fi- 

To summarize, the model has a parameter fs with dimensions of a frequency, which can 
have any value in the range < < /i, with /i ~ 1 — 100 GHz, and a dimensionless 
parameter ^ > (or equivalently /? with 2/i = |2/? — 3|). 

Performing the matching between the three phases one finally gets the spectrum. The 
exact form is 



^^gw(/) = a{p) 



(2) f af 



(27r/,)4 (f_^\ 
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(228) 



where a{fi) is a number, which depends on /u. Expanding this expression for small values 
oi f /fs one gets a result of the form 



^gw(/) 



(27r/,)4 



fs . 



ai + 



0/ 

^2 log 77^ + ^3 



(229) 



where ai, 02, 03 depend on /x. 

In the opposite limit f ^ fs (but still / smaller than the cutoff /i), one has instead 



f^gw(/) 



HiM^,^ yfs 



2^+1 



(230) 
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Figure 14: The spectrum of GWs due to amplification of vacuum fluctuations in string 
cosmology, compared with the COBE bound (dot-dashed), pulsar bound (asterisk) and 
nucleosynthesis (dotted), for a choice of the parameters fs = 10 Hz, /i = 1.5 (from ref. [l63|). 
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Figure 15: The GW spectrum, enlarging the region between the Hz and the kHz, for 
fs = 100 Hz, fi = 1.4. The dashed and dot-dashed lines are the low- and high-frequency 
limits (from ref. [|6^). 
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Both limits reproduces the frequency dependence first found in ref. |12(: , |5^. It is important 
to stress that in the high frequency limit the unknown parameter fg cancels. 

This spectrum has quite interesting properties, depending on the two free parameters 
fs and /i. First of all, it goes like at low frequencies. This result is a consequence of 
the behaviour of the cosmological model in the super inflationary pre-big-bang phase, where 
everything is under good theoretical control, so this is a rather general prediction of the 
model. This means that /iQOgw(/) can evade the COBE bound very easily, and yet be large 
in the range accessible to interferometers. This is very different from slow-roll inflation, 
where instead the spectrum is decreasing. 

If however the spectrum would behave as for all frequencies, up to the cutoff /i, 
it would be totally neglegible at, say, 1 kHz, since the value of /iQOgw(/) at the cutoff 
frequency /i is bounded by the nucleosynthesis limit discussed in sect. |7[l|. However, above 
the unknown frequency fg, the spectrum predicted by this model changes and behaves 
as Z^"^'^; this reflects the fact that the Universe changed regime and entered a phase of 
De Sitter inflation with linearly growing dilaton. The parameter measures the growth of 
the dilaton in units of the Hubble constant Hg during the string phase, and for an almost 
constant dilaton, 3 — 2// ~ 0. 

This means that, in some range for the parameters fs and fj,, it is possible to have the 
optimal situation from the point of view of interferometers: a spectrum that grows at very 
low frequencies, so that it can evade the COBE bound, and then flattens, at the maximum 
level allowed by the nucleosynthesis limit. 

The resulting spectrum is shown in fig. ffi, compared to the COBE, pulsar and nucle- 
osynthesis bounds discussed in sect. |ll|. Fig. 1^ enlarges the frequency range relevant for 
ground-based interferometers, displaying some features of the spectrum which might be a 
signature of this background. 

Finally, it is clear that in this model GWs cannot be responsible for the COBE anisotropies, 
since the ^ behaviour makes /igfigw neglegible at very low frequencies. In the pre-big- 
bang model, a possible seed for the CMBR anisotropies can be provided by pseudo-scalar 



fields |85, |3, pQ, 64, ||, 196 1. 



10 Other production mechanisms 
10.1 Phase transitions 

In the history of the Universe are expected a number of phase transition. In particular, 
the QCD phase transition takes place at T* ~ 150 MeV. Above this temperature quarks 
and gluons are deconfined and form a quark-gluon plasma, while below they are in the 
confined phase. Recent lattice calculations with three quark flavors and with the physical 
mass for the strange quark, using Wilson fermions, suggest that the phase transition is 



first order [155], but the issue is not yet settled (see also |22g| ] for a detailed discussion of 
cosmological consequencies of the QCD phase transition). Around T^, ~ 100 GeV we expect 
the electroweak phase transition, when SU{2) x U{1) breaks to U{l)em- Further phase 
transitions could occur even earlier, at the grand unified scale. These phase transitions 
are dramatic events in the history of the Universe, and are good places to look for GW 
production. 



72 



Furthermore, phase transition will in general modify previously produced spectra [234| 



This follows from the fact that during the phase transition one can have a large drop in 
the number gs of effective species (e.g. by a factor ~ 3 during the QCD phase transition). 
Since the entropy S ~ gsa^T^ is constant, this results in a smaller value of the growth rate 
of and therefore in a different rate at which modes cross into the horizon. For the QCD 
transition this results in a step in the spectrum at a frequency / ~ 10^'' Hz (determined by 
the value of the Hubble parmeter at the transition), so that the spectrum above 10"'' Hz is 
diminuished by approximately 30% [|234j| . Similar drops happen at e^e~ annihilation, with 
a step at / ~ lO"^'^ Hz, and at the GUT phase transition, with a step in the GHz region. 

10.1.1 Bubble collisions 

In first order phase transitions the Universe finds itself in a metastable state. True vacuum 
bubbles are then nucleated via quantum tunneling. In strongly first order phase transition 
the subsequent bubble dynamics is relatively simple: once the bubble are nucleated, if they 
are smaller than a critical size their volume energy cannot overcome the shrinking effect 
of the surface tension, and they disappear. However, as the temperature drops below the 
critical temperature Tc, it becomes possible to nucleate bubble that are larger than this 
critical size. These 'critical bubbles' start expanding, until their wall move at a speed close 
to the speed of light. The energy gained in the transition from the metastable state to 
the ground state is transferred to the kinetic energy of the bubble wall. As the bubble 
expands, more and more regions of space convert to the ground state, and the wall becomes 
more and more energetic. At the same time, it also becomes thinner, and therefore the 
energy density stored in the wall increases very fast. As long as we have a single spherical 
bubble, this large energy of course cannot be converted in GWs. But when two bubbles 
collide we have the condition for the liberation of a large amount of energy into GWs. In 



particular, there are two possible 'combustion' modes for the two bubbles: detonation [237|, 
that basically takes place when the boundaries propagate faster than the speed of sound, 
and deflagration, when instead they move slower. In the first case there is a large production 



of GWs III, leg, p 



The characteristic frequency, as seen today, of these GWs is given by eq. (157), with 



T^, equal to the temperature of the phase transition. First of all, we must estimate the 
parameter e that enters eq. ( |157| ). Plausible values of e have been discussed in detail by 
Hogan |148] and by Witten (ref. [^60| , app. B). In the case of the QCD and electroweak 



transitions e ~ 1 is excluded because otherwise in the bubble collisions there would be 
overproduction of primordial black holes. Production of primordial black holes is severely 
constrained because the energy density of primordial black holes scales like l/a? in the RD 
phase, and therefore they would come to dominate the energy density of the Universe. In 
particular, the Universe would not be radiation dominated during nucleosynthesis. Actually, 
primordial black holes can evaporate through emission of Hawking radiation. However, in 
order to evaporate before nucleosynthesis, the scale at which the phase transition takes 
place must be higher than 10^^ GeV p52| ]. Therefore e could be of order one in phase 
transitions which take place at these temperatures, but not in the QCD or electroweak 
phase transition. 

To estimate e it is necessary to make assumptions about how the phase transition is 
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nucleated. If it is nucleated by thermal fluctuations, ref. |14S] suggests an upper bound 
e<10^^, but recent lattice simulations (restricted however to pure QCD without fermions) 
give a much smaller value e ~ 10~^, which would make the GWs from the QCD transition 
neglegible, as will be clear below. For nucleation of bubbles via quantum tunneling a 



detailed analysis has been done in ref. 250 |. The relevant parameter is the nucleation rate; 
the phase transition can be modelled assuming an exponential bubble nucleation rate per 
unit volume, 

r = Toe^* . (231) 

The parameter /3 sets the scale for the typical time variation, and therefore the frequency 
at time of production, /^,, is determined by /?. In fact, ref. [|170(| founds that the spectrum 
peaks at 27r/* ~ 2/3. Then 

~e^K^^*- ^^^^^ 
Of course, /? is a model-dependent quantity, but for a wide class of models one typically 
finds e ~ IQ-^ - IQ-^ p8|, pOl, f7l|, 0§. 



If however the transition is nucleated by impurities (which is most often the case, except 
in very pure and homogeneous samples) the issue is much more complicated. For instance, 
the impurities could be given by turbolent motion of the cosmic fluid generated prior to the 
QCD epoch; they would depend on the detailed spectrum of the fluid motion, or in general 
on the characteristic distance between impurities, and in this case e is basically impossible 



to estimate, even as an order of magnitude | 260 |. 



Coming back to the simplest case of nucleation by quantum or thermal fluctuations, for 



the electroweak phase transition, ref. [170| finds 



^ ~ 1.3 X 10^3 (233) 

and therefore 

/o ~ 4.1 X lO^^Hz. (234) 

This phase transition is therefore very relevant for LISA. With the same value for e, the 
QCD phase transition would instead peak around /o ~ 4 x 10^^ Hz. A smaller value of e is 
required to bring it into the frequency window of LISA. 

A first estimate of the value of /igOgw, that as we will see really holds only for strongly 
first order phase transitions, can be given as follows |252| ]. The typical wavelength A,, 
produced in the collision of two bubbles will be of the order of the radius of the bubbles 
when they collide, which in turn is a fraction of the horizon scale H~^. The energy liberated 
in GWs in the collision of two bubbles is of order 

Egw ~ , (235) 

where Mb ~ /OvacA^ is the energy of a typical bubble and /Ovac the difference in energy 
density from the metastable state to the ground state. The fraction of the false vacuum 
energy that goes into GWs instead of going into is therefore Eq^^ /Mb- Since 

E'gW GMb 2 /OQftA 

-irj 7 Gp^^cK (236) 

Mb a* 
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Gpvac, we get, at time of production 

Pgw ~ >^Ih^Pj = f^P-y ■ (237) 

Therefore pgw is related to p^y, and the basic reason is that there is essentially only one 
dimensionful parameter that enters in the estimate. If we neglect the variation in g{T) 
from time of production to the present value, Pgw/P'y is almost constant, and since today 
hoP-y/Pc — 2.5 X 10~^, this quick order of magnitude estimate gives /io^gw(/o) ~ 10~^e^. 



However, further suppression factors appears in the accurate computation of ref. [|170| , 
which gives 

Here a is the ratio of vacuum energy to thermal energy in the symmetric phase. It charac- 
terizes the strength of the phase transition; for a — > we have a very weak first order and 
for a — > oo a very strong first order phase transition. For the electroweak phase transition, 
ref. |17|] gives 

a = 1.4 X 10"^ . (239) 

K is an efficiency factor quantifying the fraction of vacuum energy which goes into bubble 
wall kinetic energy, and in the electroweak case is estimated to be ac ~ 7.8 x 10~^. Finally, v 
is the velocity of the bubble wall, and is approximately equal to the speed of sound during 
RD, V ~ l/\/3. Altogether, these suppression factors give a disappointing 

/igjlgw ~ 10^22 . (240) 

These suppression factors are a consequence of the fact that the electroweak phase transition 
is very weakly first order, if first order at all. It is clear that strongly first order phase 
transitions are needed for a detectable signal. 

Quite interestingly, the requirement of strongly first order phase transition is the same 
that it is necessary to make possible electroweak baryogenesis (see e.g. ref. |220| ] for a recent 
review). In particular, the ways in which baryons can be produced have been separated 
into 'local baryogenesis', when both baryon number violating and CP violating processes 
occur near the bubble wall, and 'nonlocal baryogenesis', when only CP violating processes 
take place near the wall. The condition for local baryogenesis to dominate is that the speed 
of the wall be greater than the speed of sound, which is the same condition for detonation 
to dominate the bubble-bubble collisions. 

Partly motivated by the importance for baryogenesis, there have been many recent 
investigations on the strength of the phase transition in extensions of the Standard Model, 



and in particular in the Minimal Super symmetric Standard Model (MSSM) [ll34| , 203 , 103, 



Both analytical pB, 97, 102, BO, 183, 107, m and lattice WTE, nm, m computations 



have shown that in the MSSM the phase transition can be strongly first order if the top 
squark is lighter then the top quark. The observation of a relic stochastic background at 
LISA could therefore give us a hint for physics beyond the Standard Model. 
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10.1.2 Turbolence 



For weakly first order phase transitions the situation is actually more complex. Bubble 
nucleation occurs both via quantum tunneling and with thermal fluctuations, and the sub- 
sequent evolution of the bubble depends on the interaction of the wall with the surrounding 
plasma; part of the energy gained in the transition from the metastable state to the ground 
state is used to heat up the plasma, and another fraction is converted into bulk motion of 
the fluid. If the Reynolds number of the Universe at the phase transition is large enough, 
then this results in the onset of turbolence in the plasma, and this is another powerful 
source of GWs. 

The estimates for the characteristic frequency and for h^Qg^ turn out to be [170| 



/o 2.6 X 10-« Hz ^ f AU ( ^) , (241) 



Hj \lGe\) UOOy 
ht^g^ ^ 10-^ [^) ' vvl (—) , (242) 



where vq is the fluid velocity on the largest length scales on which the turbolence is being 
driven, and v the velocity of the bubble wall. These estimates indicate that fully developed 
turbolence can be comparable, or even more potent than bubble collisions in generating 
GWs. 

10.1.3 Scalar field relaxation 

Consider a global phase transition in the early Universe, associated with some scalar field 
(j) that, below a critical temperature, gets a vacuum expectation value {(j)). For causality 
reasons, in an expanding Universe the scalar field cannot have a correlation length larger 
than the horizon size. So, even if the configuration which minimizes the energy is a constant 
field, the field will be constant only over a horizon distance. During the RD phase the 
horizon expands, and the field will relax to a spatially uniform configuration within the 
new horizon distance. This relaxation process will in general produce gravitational waves, 
since there is no reason to expect (T^i/) for the classical field just entering the horizon to 



be spherically symmetric [173]. A simple estimate of the spectrum produced makes use of 
the fact that the characteristic scale of spatial or temporal variation of the field is given by 
the inverse of the Hubble constant, H~^^ at the value of time under consideration. Then, 
the energy density of the field is p ~ (dcp)'^ ~ The total energy in a Hubble volume 

is pH~^ ~ {(j))'^H~^ and the corresponding quadrupole moment is Q ~ {pH~^)H~^ ~ 
{cj))'^H^^, times a constant smaller than one, which measures the non-sphericity of the 



configuration. The energy liberated in GWs in a horizon time is then given by |]17 



A^; ~ H-^ X Luminosity ~ R-^G { -j-f ] ~ G{(p)^H~^ , (243) 



and the energy density of GW produced is 
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This is the energy density at a frequency / ~ -ff, at a value of time t = which is 

the time when the mode with frequency / enters the horizon. The asterisk in {dp/d\og /)* 
reminds that this quantity is evaluated at horizon crossing, therefore at different values of 
time for different frequencies. However, in the RD phase, the energy density scales like 
p ~ l/a^(t) ~ and H ~ so that p scales as H^. Therefore, eq. ( |244| ) means that, if 
we evaluate dp/dlogf at the same value of time, for all frequencies which at this value of 
time are inside the horizon, the spectrum is flat, /igOgw(/) ~ const. This is a rather nice 
example of how a flat spectrum can follow simply from the requirement of causality and 



simple dimensional estimates. The simplest way to redshift eq. ( 244 ) to the present time is 
to use = {8Tr/3)Gpy, so that eq. ( p44D can be written as 



c^log//* 3 



G\<P)Hp,)*, (245) 



where both sides are evaluated at horizon crossing. Then pg^ redshifts as p-y adia- 

batically, eq. (153), and p^a.d = 9{T)P'y, so that (r2gw)o = {do/ g*Y^^{^sw)*) and therefore 



eq. ( 244 ) gives, for the energy density today, 

1/3 



JO .5xlO-f#Vf™V'. (246) 



gw 

A quite similar mechanism has been studied in ref. [|149(| , considering a theory with a 
multicomponent scalar field cj) and an effective potential V{(j)) with degenerate minima 
at a set of points with \(j)\ = (pQ ^ 1. The direction in field space with dV/d(j) define 
massless Goldstone bosons. Cooling down from a state with higher symmetry generates 
spatial gradients in (j) and therefore excites these Goldstone bosons. The modes that enter 
the horizon then produce GWs with the same mechanism described above, and again one 
finds eq. ( pie] ). 

Clearly, the effect is significant only if {(j)) is not too far from the Planck mass. Since 
the mechanism is completely general, it is well conceivable that it might take place in string 
theory, which has a multitude of scalar fields, like the dilaton and the moduli of the com- 
pactification. At tree level these fields correspond to fiat directions, but non-perturbative 
eff'ects will lift the degeneracy; it is then very natural to expect that at the same time they 
produce a number of degenerate or almost degenerate vacua. In this case this effect could 
give an extremely interesting signal both for LISA and for ground based interferometers. 

Of course, if we want to avoid the COBE bound, it is necessary to suppress the spectrum 



at low frequencies. As suggested in ref. |14£], this might happen if the Goldstone modes are 
coupled to fields that strongly damp their oscillations after a certain epoch. For example 
there could be a second phase transition that removes the degeneracy between the minima 
of the potential; the field would then relax everywhere to a single minimum, suppressing 
further GW production. 

Finally, another related production mechanism is the collapse of solitons or of unstable 
domain walls [135, p.36| ]. In this case the spectrum is quite peaked, at a frequency determined 



by the expectation value of the scalar field and for of the order of the electroweak 
scale the signal would be relevant for LISA. 
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10.2 Cosmic strings 

Cosmic strings are topological defect that might have formed during phase transitions in 



the early Universe (see e.g. |254, |255(] and references therein). Interest in cosmic strings 
has also been fueled by the fact that they were considered a candidate for seeding structure 
formations. The most recent comparison with CMB anisotropics show that in models with 
vanishing cosmological constant they do not reproduce the dependence of the anisotropy 
from the multipole moment |214| , |6|. However, with a non- vanishing cosmological constant 
they can still be a viable option p^, IM] . 



The spectrum of GWs produced by cosmic strings is discussed in detail in refs. [^54 , 



255| , H, |, 1^, ||. Gravitational waves are produced by the relativistic oscillations of a 



cosmic string, in a cosmic string network. The dynamics of this string network has been 
studied in detail, and it has been found that it obeys a scaling property, so that the only 
relevant scale is given by the Hubble length. Small loops oscillate producing GWs and 
disappear, while more small loops are chooped of very long strings, replacing the loops that 
disappear, and the network has a self-similar configuration, with loops of all lengthscales. 
Since the wavelength of the GWs emitted is fixed by the length of the loop, the spectrum 
of the cosmic string network extends across a very large frequency band. 

The spectrum has two main features: an almost flat region that extends from / ~ 10~^ 
Hz up to / ~ 10^'' Hz, and a peak in the region / ~ 10~^^ Hz. 

At LIGO/VIRGO frequencies we are in the flat part of the spectrum, and the typical 
estimate of the intensity is of order 

/i^figw ~ 10"^ - 10"^ . (247) 



(Similar results are obtained from hybrid topological defects |192(] ). The scale for these 
numbers is given first of all by the combination (G/u)^, where is the mass per unit length 
of the string. It comes from the fact that a loop radiates with a power P ~ ^Gfi'^, (7 is a 
dimensionless constant) while, evaluating h^Qg^, another factor of G comes from 1/pc- For 
strings created at the GUT scale G// ~ 10~^. Then, there are also large numbers related to 
the number of strings per horizon volume, which is thought to be of order 50, to the fact 
that also 7 ~ 50, and to the size of the loops at formation time ||8|, that finally give a value 
/igJlgw ~ 10-8 - 10-"^. 

For the cosmic string spectrum the most relevant bound comes from the pulsars, dis- 
cussed in sect. Pulsar timing implies a bound on the mass per unit length of the 



string [g3| 



^ < (5.4±1.1) X 10"^ (248) 
For comparison, the value of obtained fitting CMBR anisotropics in an open Universe is 



{Gfi)/c ^ 1.7 X 10 for i7o = 0.2 [27|. The improvement in pulsar timing expected in the 



next few years, see sect. 11, will therefore be crucial for testing the cosmic string scenario. 



10.3 Production during reheating 

If we have an inflationary phase below the Planck scale, there is no gravitational analog of 
the 2.7K radiation. At the Planck scale photons and gravitons can be produced simply by 
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thermal collisions, with probably similar production rates. However, the particles produced 
in this way at the Planck (or string) scale are exponentially diluted by the subsequent 
inflationary phase, and the photons that we see today in the cosmic microwave background 
radiation (CMBR) have been produced during the reheating era which terminated the 
inflationary phase. 

Since the reheating temperature Ti-h < Minfl <^ Mpi, thermal collisions are by now 
unable to produce a substantial amount of gravitons. However, in this case relic GWs 
can be produced through some non-equilibrium phenomenon connected with the reheating 
process. For a process that takes place at the reheating temperature Tj-h the estimate of 
the characteristic frequency is given by eq. ( |157| ) with = T^h. In principle the reheating 
temperature can be between a minimum value at the TeV scale (since this is the last 
chance for baryogenesis, via the anomaly in the electroweak theory) and a maximum value 
Trh ~ Minfl<3 X 10^^ GeV. In supersymmetric theories, the gravitino problem give a further 
constraint T^t^W^ GeV Il6|,[l6|]. More precisely, ref. |6|] gives a bound on Tj-h ~ 10^ — 10^ 
GeV for a gravitino mass 100 GeV-1 TeV, and a bound lO" - lO^^ Qgy independent of 
the gravitino mass. 

So, phenomena occurring at reheating can manifest themselves with cutoff frequencies 
between 10~^ and 10^ Hz, if we set e = 1 in eq. ( |157| ). In typical non-super symmetric models 
the reheating temperature is large, say 10^^ GeV, corresponding to /o ~ (l/e)10^ Hz, and 
in this case interferometers or resonant masses could only look for their low-frequency 
tails. One should however keep in mind the possibility that bubble nucleation occurs before 
the end of inflation, and then the cutoff frequency would be redshifted by the subsequent 
inflationary evolution toward lower values, possibly within the VIRGO/LIGO frequency 
range [ pO[ . 

If instead the reheating temperature is ~ 10^ — 10^ GeV, non-equilibrium phenomena 
taking place during reheating would give a signal just in the LIGO/VIRGO frequency band. 
A reheating temperature T ~ TeV would instead be relevant for LISA. 

A mechanism for GW generation during reheating is bubbles collision, in the case when 



inflation terminates with a first order phase transition [ 252 , |17Cf| . We have already discussed 



it in sect. 10.1.1 



Another possibility is that reheating occurs through the decay of the inflaton field. In 
this case, in a very general class of models, there is first an explosive stage called preheating, 
where the inflaton field decays through a non-perturbative process due to parametric reso- 
nance 1 167], and at this stage there are mechanisms that can produce GWs, see [166, 32 1. 
We believe that in these cases a reliable estimate of the characteristic frequency is really 
very difficult to obtain, since the relevant parameter, which is the value of the Hubble con- 
stant at time of production, depends on the complicated dynamics of preheating and on the 
specific infiationary model considered. Ref. |166f| examines models where the characteristic 
frequency turns out to be in the region between a few tens and a few hundreds kHz. 



11 Stochastic backgrounds of astrophysical origin 

The emission of GWs from a large number of unresolved astrophysical sources can create a 
stochastic background of GWs. This background would give very interesting informations 
on the state of the Universe at redshifts z ~ 2 — 5. It will provide a probe of star formation 
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rates, supernova rates, branching ratios between black hole and neutron star formations, 
mass distribution of black holes births, angular momentum distributions and black hole 
growth mechanisms |^^. At the same time, however, from the point of view of cosmological 
backgrounds produced in the primordial Universe, the astrophysical background is a 'noise', 
which can mask the relic cosmological signal. 

In this section we investigate a number of astrophysical mechanisms that can produce 
a stochastic background of GWs. As for a general orientation into the typical intensity 
and frequency, we should note that GWs of astrophysical origin of course are not subject 
to the nucleosynthesis bound, since they were created much later, so that one of our main 
'benchmarks' disappears in this case. The estimates of the typical values of /ig^^gw is strictly 
dependent on the mechanisms, and we will see examples below. 

Concerning the frequency, a first observation is that there is a maximum frequency at 
which astrophysical sources can radiate. This comes from the fact that a source of mass M, 
even if very compact, will be at least as large as its gravitational radius 2GM, the bound 
being saturated by black holes. Even if its surface were rotating at the speed of light, its 
rotation period would be at least 4ttGM, and the source cannot emit waves with a period 
much shorter than that. Therefore we have a maximum frequency [ p45| , 

To emit near this maximum frequency an object must presumably have a mass of the 
order of the Chandrasekhar limit ~ 1.2Mq, which gives a maximum frequency of order 



10 kHz p45|| , and this limit can be saturated only by very compact objects (see ref. [ 110 | 
for a recent review of GWs emitted in the gravitational collapse to black holes, with typical 
frequencies kHz). The same numbers, apart from factors of order one, can be obtained 
using the fact that for a self-gravitating Newtonian system with density p, radius R and 



mass M = p(4/3)7ri?^ there is a natural dynamical frequency 123 1[| 



A„.-L„G,)V^.(^)"\ P50) 

With R > 2GM we recover the same order of magnitude estimate apart from a factor 
(3/8)1/2 ^ 0.6. This is an interesting result, because it shows that the natural frequency 
domains of cosmological and astrophysical sources can be very different. In particular, a 
GW signal detected above, say, 10 kHz, would be unambiguosly of cosmological origin. 
However, as we have seen, it is not easy to build sensitive detectors operating at such high 
frequencies. 

In discussing a stochastic background of astrophysical origin, apart from the character- 
istic frequency and characteristic intensity, there is a third important parameter, which, 
for burst sources, is the duty cycle D. This is given by the ratio of the typical duration 
of the signal (e.g. 1 ms for supernovae) to the average distance between successive bursts. 
If L> <C 1, the signal is not stochastic. For intermediate values, say D ~ 0.1, we have a 
so called 'popcorn' noise, or shot noise. As D — > 1, we get a continuous, stochastic back- 
ground. The value of D is important because, if D <^ 1, we can distinguish a signal from 
the very early Universe form an astrophysical signal, even when the latter corresponds to 
higher peak amplitude. 
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11.1 Supernovae collapse to black-hole 

Sufficiently massive stars, at the final stage of their evolution, collapse to form a black hole or 
a neutron star. In this explosive supernova event GWs are liberated. Depending on the rate 
of supernova events, this can give origin to a stochastic background. The typical timescale 
of supernovae is the millisecond, and therefore the typical frequency of the GWs produced 
is around the kHz. These processes are therefore relevant for ground based experiments. 



Earlier estimates |47, 45] suggested a duty cycle D ~ 0.3 and a value /igilgw possibly 
even of order 10~^ at the kHz, but acknowledged the need for better collapse models and 
duty cycle informations. 



Recently, a careful study has been done in ref. |108|. The crucial point is the deter- 
mination of the star formation rate. However, in the last few years the understanding of 
the origin and evolution of galaxies has greatly improved, thanks to the spectacular data 
from the Hubble Space Telescope, Keck, and other large telescopes, that have allowed to 
investigate the Universe up to redshifts z ~ 4 — 5. This makes possible a determination 
of the star formation rate density evolution, based on observations. Denote by the 
comoving star formation rate density, i.e. the mass of gas that goes into stars, per unit time 
and comoving volume element, at redshift z. Then the rate of core-collapse supernovae 
event Rs^{z), i.e. the number of events per unit time in a comoving volume, at redshift z, 
is |0| 

Rsn{z) = dz' ^ ^M)dM, (251) 

Jo dz' Z' + L JMp 

where $(M) is the Salpeter initial mass function and dV/dz the comoving volume element. 
The lower limit Mp depends on the specific nature of supernovae considered: the smallest 
progenitor mass which is expected to lead, after collapse, to black hole is in the range 
(18 - 30) M0, and the reference value 25M0 has been used in |lO|], while M„ = 125M0. 



With these values, adopting a value /iq = 0.5 for the Hubble constant and assuming a flat 
cosmology, VLq = 1, with vanishing cosmological constant, the total number of supernovae 



explosions per unit time leading to black hole formation is |108| 

i?BH = 4.74 events/s . (252) 
The duty cycle due to all sources up to redshift z is defined as 

Diz)= r dRBuATGwil + z), (253) 
Jo 

where Atgw is the average time duration of single bursts at the emission, and is of order 



1ms. Then the total duty cycle turns out to be [108] 



D = 1.57 X 10^2 . (254) 

This implies that this GW background is not stochastic, but rather a sequence of bursts, 
with typical duration 1ms and a much larger typical separation ~ 0.2 s. This is a good new 
from the perspective of primordial Universe cosmology, since this astrophysical background 
can be distinguished from relic GWs produced at much higher redshifts. 



The value of /ig^gw has been computed in ref. ]105] and is of order 10 peaked 
around 2 kHz. It is shown in fig. |l^ (provided by the authors of ref. JlOS], and including a 
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Figure 16: rjgw against frequency, for GWs emitted in the supernovae collapse to black 



holes, from ref. |10S 



recent update of the star formation rate at large z) for three different values of the rotation 
parameter a = Jc/{GM^^j.^), where J is the angular momentum. A black hole forms only if 
a is smaller than a critical value estimated in the range acrit ~ 0.8 — 1.2; otherwise rotational 
energy dominates, the star bounces, and no collapse occurs. For comparison with the other 
plots of /igfigw shown in this report, note that fig. 16 shows Jlgw rather than h^^lgw, and 



that ref. |108| uses /iq = 0.5, so the value of /igilgw is obtained from fig. ^ dividing by 4. 
Note also that both axes are on a linear scale, rather than a log-log scale, as the other plots 
that we show. 

From the figure, we see that this background is out of reach for first generation experi- 
ments, while it could be detectable with advanced ground based interferometers. 

It is also important to observe that the resulting GW signal is insensitive to the uncer- 
tainties that are present in the star formation rate in the high redshift region z ~ 4 — 5, 
since the result turns out to be dominated by low-to-intermediate redshift sources. 

11.2 GWs from hydrodynamic waves in rotating neutron stars 

In rapidly rotating neutron stars there is a whole class of instabilities driven by the emis- 
sion of GWs, called CFS instabilities after Chandrasekhar ffsll, who discovered them, and 



Friedman and Schutz |122], who showed that they are generic to rotating neutron stars. 
Basically, they form when there is a mode that is forward-going, as seen from a distant 
observer, but backward going with respect to the rotation of the star. In this case, when 
this mode radiates away angular momentum, the star can find a rotation state with lower 
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Figure 17: /igOgw against frequency, for GWs produced by the r-mode of rotating neutron 



stars (from ref. |109|). 



energy, than fueling further growth of the mode and therefore an instability occurs. Viscous 
forces, instead, tend to damp this instability. 

Recent studies of rapidly rotating relativistic stars have revealed the existence of a 
particularly intersting class of modes (r-modes) that are unstable due to the emission of 
GWs [14, 121, |18C| ]. As the star spins down, for these modes an energy of order of 1% of a 
solar mass is emitted in GWs, making the process very interesting for GW detection [ ^07| , 



15|. Ref. |207| give a rough estimate of the spectrum assuming a comoving number density 
of neutron star births constant in the range < z < 4, and zero at higher z. A more detailed 



analysis has been done in ref. [109, 230 [, using again the star formation rate determined by 
observation, mentioned in the previous section. The analysis is performed for three different 
cosmological models, a flat Universe with vanishing cosmological constant (f^Af = Ij = 0) 
and ho = 0.5, a flat, low-density model {^m = 0.3, I^a = 0.7, ho = 0.6), and an open model 
with Qm = 0.4, J7a = 0, /iq = 0.6. Similarly to the computation of R^s in the previous 
section, one can now compute the number of neutron stars formed per unit time within 
a comoving volume out at redshift z. Depending on the cosmological model, and on the 
upper cutoff in the progenitor mass, it turns out that 



i?NS — (15 — 30) events/s . 



(255) 



Comparing with the collapse to black hole, with the same parameters, one finds that i?Ns ~ 

4i2BH. 

The mechanism of GW emission is discussed in refs. [|180| , ^07| , 16 [. The neutron 
star is modeled as a fluid with a polytropic equation of state p = kp^, with k chosen so 
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that the mass of the neutron star is IAMq and the radius is R = 12.53 km. The r-mode 
excitation starts as a small perturbation of the velocity field. As long as the amplitude 
a of the r-mode us small, the effect of the instability is such that a grows exponentially 
while the angular velocity Q of the neutron star is nearly constant. After approximately 
500s, a becomes of order one, and we reach a regime where non-linear hydrodynamic effects 
are important. The non-linear effects saturate the growth of a, which stays approximately 
constant. This phase lasts approximately for a time tns ~ 1 yr, during which the star 
loses approximately 2/3 of its rotational energy by emission of GWs. Finally, the angular 
velocity becomes sufficiently low so that the r-mode ceases to be unstable. 

Since GWs are emitted continuously over approximately lyr, the duty cycle for this 



process is [ll09|] 



D= I dRmfm{l + z) ~ 10^ (256) 







and therefore we definitely have a stochastic background. 

The minimum and maximum values of the frequency radiated are estimated as fol- 
lows [207|. It is known that the maximum possible angular velocity of a rotating star 
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1/2 

(257) 



where M, R are the mass and radius of the corresponding non-rotating star. Assuming that 
most neutron stars are born with close to Oi^, one can write the initial rotational energy 
of the star. Using the fact that about 2/3 of it is radiated into GWs, one finds 

/max ^^^K. (258) 

ovr 

The lower cutoff in / is associated to the fact that the r-mode instability becomes ineffective 
below some critical rotational frequency, and it can be estimated that 

/^ax=^ 120 Hz. (259) 



The resulting spectrum of GWs is shown in fig. |T7|, (provided by the authors of ref. [ p.09| , 
and including a recent update of the star formation rate at large z and other technical 
improvements) for the three different cosmological models defined above. We see that the 
three curves are very similar, and therefore the result is quite solid against changes in the 



cosmological model, and [109 



/lo^gw =^ 8 X 10"^ (260) 

in the frequency range (500-700) Hz. This value is quite interesting, but from the discussion 
in sect. |6|it is clear that it is still below the sensitivities that can be obtained with correlations 
of first generation interferometers. It will be instead quite accessible to second generation 
experiments. 

11.3 GWs from mass multipoles of rotating neutron stars 

In the previous subsection we examined the case where GWs are emitted by rotating neutron 
stars through their coupling to the current multipoles associated to hydrodynamic waves 
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on the star surface. We now discuss the emission by the simplest mechanism, which is the 
couphng of the GWs to the multipoles of the mass distribution of a rotating neutron star. 

The stochastic background from mass multipoles of rotating neutron stars has been dis- 
cussed in ref. |217]. The main uncertainty comes from the estimate of the typical ellipticity 



e of the neutron star, which measures its deviation from sphericity. An upper bound on e 
can be obtained assuming that the observed slowing down of the period of known pulsars 
is interely due to the emission of gravitational radiation. This is almost certainly a gross 
overestimate, since most of the spin down is probably due to electromagnetic losses, at least 
for Crab-like pulsars. 

The typical frequency for these GWs can be in the range of ground-based interferometers. 



With realistic estimates for e, ref. [217] gives, at / = 100 Hz, a value of /ic(/) 5 x 10 
that, using eq. (|^), corresponds to 

/igl7gw(100Hz) ~ 10"^^ . (261) 

This is very far from the sensitivity of even the advanced experiments. An absolute upper 
bound can be obtained assuming that the spin down is due only to gravitational losses, and 
this gives /i^J7gw(100Hz) ~ 10-^ but again this value is probably a gross overestimate. 

Techniques for the detection of this background with a single interferometer using the 
fact that it is not isotropic and exploiting the sideral modulation of the signal have been 



discussed in ref. 1133 



11.4 Unresolved galactic and extragalactic binaries 

These backgrounds will be relevant in the LISA frequency band. For frequencies below a 
few mHz, one expects a stochastic background due to a large number of galactic white dwarf 



binaries []T05|, l245|, pT], ^1% m 



In sections ( 11.2 . 11.1 ), in order to determine whether the superposition of burst signals 
from supernovae formed a stochastic background, we used the duty cycle. For a continuous 
signal as that produced by binaries, instead, whether we have a stochastic background or 
not depends on our frequency resolution. After a time T of observation, we can resolve 
A/ = 1/T; for T = 1 yr, A/ ~ 3 x 10~® Hz, and the number of resolvable frequencies in 
the LISA frequency band is of order (lHz)/(3 x 10~^ Hz) = 3 x 10^. It is estimated that 
most frequency bins below a critical frequency of the order of 1 mHz will contains signal 
from more than one galactic binary, and will therefore form a confusion- limited background. 
Above this critical frequency, instead, individual signals will be resolved. 

In particular, it is expected a contribution from compact white dwarf binaries. The 
computation of the intensity of this background depends on the rate of white dwarf mergers, 
which is uncertain. A possible estimate is shown in fig. ^ ||3^. We see that, below a few mHz, 
it would cover a stochastic background of cosmological origin at the level h^Qgw ~ 10^^^. 

It should be observed that, even if an astrophysical background is present, and masks 
a relic background, not all hopes of observing the cosmological signal are lost. If we un- 
derstand well enough the astrophysical background, we can subtract it, and the relic back- 
ground would still be observable if it is much larger than the uncertainty that we have on 
the subtraction of the astrophysical background. In fact, LISA should be able to subtract 
the background due to white dwarf binaries, since there is a large number of binaries close 
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enough to be individually resolvable jS^]- This should allow to predict with some accu- 
racy the space density of white dwarf binaries in other parts of the Galaxy, and therefore 
to compute the stochastic background that they produce. Furthermore, any background of 
galactic origin is likely to be concentrated near the galactic plane, and this is another handle 
for its identification and subtraction. The situation is more uncertain for the contribution 
of extragalactic binaries, which again can be relevant at LISA frequencies. The uncertainty 
in the merging rate is such that it cannot be predicted reliably, but it is believed to be lower 
than the galactic background |217|. In this case the only handle for the subtraction would 
be the form of the spectrum. In fact, even if the strength is quite uncertain, the form of 
the spectrum may be known quite well p6[| . 

Finally, another possible astrophysical background can be expected if MACHO's are 
black holes of mass ~ 0.5Mq |204|. In this case we can expect that a large number of them 
will be in binary systems, and they could create a stochastic background relevant both for 
LISA and for ground-based interferometers |154]. 



12 Conclusions 

Present GW experiments have not been designed especially for the detection of GW back- 
grounds of cosmological origin. Nevertheless, there are chances that in their frequency 
window there might be a cosmological signal. The most naive estimate of the frequency 
range for signals from the very early Universe singles out the GHz region, very far from the 
region accessible to ground based interferometers, / < a few kHz, or to resonant masses. To 
have a signal in the accessible region, one of these two conditons should be met: either we 
find a spectrum with a long low-frequency tail, that extends from the GHz down to the kHz 
region, or we have some explosive production mechanism much below the Planck scale. As 
we have discussed, both situations seem to be not at all unusual, at least in the examples 
that have been worked out to date. The crucial point is the value of the intensity of the 
background. 

With a very optimistic attitude, one could hope for a signal, with the maximum intensity 
compatible with the nucleosynthesis bound, /ig^^gw ~ a few xlO^^ (or even 10~^, stretching 
all parameters to the maximum limit). Such an option is not excluded, and the fact that 
such a background is not predicted by the mechanisms that have been investigated to date 
is probably not a very strong objection, given our theoretical ignorance of physics at the 
Planck scale and the rate at which new production mechanisms have been proposed in recent 
years, see the reference list. However, with more realistic estimates, on general grounds it 
appears difficult to predict a background that, either in the region between a few Hz and a 
few kHz, relevant for ground based interferometers, or in the region (10~^ — 1) Hz relevant to 
LISA, exceeds the level hQ^lg^ ~ a few xlO~^, independently of the production mechanism. 
This should be considered the minimum detection level where a significant search can start. 
Such a level is beyond the sensitivity of first generation ground based detectors, but it is 
accessible both to advanced interferometers and to the space interferometer LISA. 

The difficulties of such a detection are clear, but the payoff of a positive result would be 
enormous, opening up a window in the Universe and in fundamental high-energy physics 



86 



that will never be reached with particle physics experiments. 
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